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Abstract 

Casimir effect in most general terms may be understood as a backreaction 
of a quantum system causing an adiabatic change of the external conditions 
under which it is placed. This paper is the second installment of a work 
scrutinizing this effect with the use of algebraic methods in quantum the- 
ory. The general scheme worked out in the first part is applied here to the 
discussion of particular models. We consider models of the quantum scalar 
field subject to external interaction with "softened" Dirichlet or Neumann 
boundary conditions on two parallel planes. We show that the case of electro- 
magnetic field with softened perfect conductor conditions on the planes may 
be reduced to the other two. The "softening" is implemented on the level of 
the dynamics, and is not imposed ad hoc, as is usual in most treatments, on 
the level of observables. We calculate formulas for the backreaction energy 
in these models. We find that the common belief that for electromagnetic 
field the backreaction force tends to the strict Casimir formula in the limit 
of "removed cutoff" is not confirmed by our strict analysis. The formula is 
model dependent and the Casimir value is merely a term in the asymptotic 
expansion of the formula in inverse powers of the distance of the planes. Typ- 
ical behaviour of the energy for large separation of the plates in the class 
of models considered is a quadratic fall-of. Depending on the details of the 
"softening" of the boundary conditions the backreaction force may become 
repulsive for large separations. 
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1 Introduction 



This is the second of the two papers devoted to the Casimir effect in which we 
develop more fully what was announced in T. In these papers we advocate the 
use of the algebraic approach to the quantum systems as the most natural setting 
for the discussion of the Casimir effect. This approach gives a clear understanding 
of the sources of the difficulties one encounters in more traditional treatments, 
and allows a mathematically rigorous analysis of the effect. In the first paper P 
this analysis was carried out on a more general level, for a wide class of quantum 
systems and external conditions. We derived a general criterion for admissible 
models and obtained formulae for the backreaction energy and generalized force. 
The reader should refer to P for the background and our statement of the problem, 
and the results mentioned here. Sections 1 - 5 of that paper form prerequisites 
for the present one, and in what follows we assume their knowledge by the reader. 
Also, we refer the reader to I for bibliography. 

The main results and outline of the paper are as follows. In Section|21we show 
how the quantum scalar field with external interaction of planar symmetry fits into 
the scheme discussed in Sections 1-3 and 1-5. More precisely, the perturbing inter- 
action is assumed to modify the z-motion only (the system is then translationally 
symmetric in the directions orthogonal to the z-axis). The necessary preliminary 
condition (ii) of Sec. 1-5 for the admissibility of the model is formulated in terms 
of the z-motion generator. In Section |21 we obtain the necessary and sufficient 
condition for the equivalence of the vacuum representation and the ground state 
representation of the field influenced by the perturbation. In Section^ this con- 
dition is then extended to be also necessary and sufficient for the existence of 
finite backreaction energy. Next, in Section Owe propose a class of models of the 
field perturbation imitating the Dirichlet or Neumann conditions on two parallel 
planes separated by an adiabatically changeable distance. The models depend on 
two functions whose role is to soften the effect of the planes for high energies of the 
particles. For a class of these functions the admissibility conditions are satisfied 
and an explicit formula for the Casimir energy is obtained. We investigate this for- 
mula in Sectional We show that the original Casimir expression constitutes a term 
in the expansion of the formula in inverse powers of the distance of the planes. 
Using a scaling property of the formula we can approximate the strict boundary 
conditions. In the case of Dirichlet conditions the Casimir energy becomes infinite, 
but the Casimir force may be interpreted to approach the Casimir value in the 
limit (although with qualifications). However, in the Neumann case both energy 
and force become meaningless in the limit. We view this as a typical situation, the 
Dirichlet case being exceptional. Section [7| treats the electromagnetic case with 
metallic boundary conditions. We show that this model may be reduced to the 
superposition of scalar Dirichlet and Neumann cases. Therefore, there is no strict 
boundary limit for this model. We find that in the class of models considered here 
typical fall-off of the backreaction force for large separation of the plates is by one 
order weaker than in the original Casimir formula, and the force may become re- 
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pulsive in the limit. Appendices contain more technical results needed in the main 
text. 

Precise formulation of the central results of the paper is given in (Mod) in 
Section |31 (Asym) in Sectional and in Sectional Also, more extensive discussion 
of their physical significance will be found in the opening parts of Sections [S] and El 
and in the closing part of Section [7| 

2 Scalar field under external conditions with pla- 
nar symmetry 

We apply here the formalism of Sections 1-3 and 1-5 to more specific models. We 
take the real Hilbert space TZ of Sec. 1-3 to be the tensor product of a space 
TZ± of the motion in the x-y plane, and of a space TZz of the motion in the z- 
direction, TZ = TZ± (g) TZz- For the complexified versions of these spaces we also 
have JC ~ IC± <S) JCz- Let a positive operator h± in IC±, with domain 'D{h±), 
describe the perpendicular motion, and a positive operator hz in K-z , with domain 
Vihz): describe the z-motion. (Both operators are assumed to commute with the 
complex conjugation.) Then the operator h in K, defined in standard way by the 
form method as 

v{h)=v{h^®id)r\V(;id®hz), h = y/{h^®idf + {;id®hzf , (2.1) 

is a densely defined, selfadjoint, positive operator. More precisely, the prescription 

(?(V',¥') = ((/ij^(g)id)?/', (ft_L«)id)v3) + {(id®hz)i^,iid®hz)'^) (2.2) 

defines a closed quadratic form on 'D{h± (g) id) n 'D{id^hz)- This form determines 
in the standard way the unique positive operator h by q{Tp,(p) — (hip, hip). This 
operator defines a quantum model, as described in Sec. T3. 

The free quantum scalar field model fits into this scheme with the choices 
/C_L = L^{M.^,dxdy), ICz = L^(R,dz), h^ = -A_l, hi = -d^ where A_l is the 
two-dimensional Laplacian. We now want to introduce external conditions which 
modify the z-motion, while leaving the transversal motion unchanged. However, if 
we leave the strict translational symmetry in the x-y plane intact, then the bound 
(T5.8), which gives the condition for the implementability of the modified dynamics 
in the original representation (technically: the implementability of an appropriate 
Bogoliubov transformation), cannot be satisfied. Both mathematical and physical 
reason is clear: this condition states that in the new state the particle number A/'a 
is finite, which cannot be satisfied due to the translational symmetry in the x-y 
plane. Nevertheless, we can apply the standard "thermodynamic limit" procedure 
for the transversal directions: we restrict the x-y motion to some compact region, 
demand finiteness of particle number A/'a and Casimir energy £a (conds. (1-5.8,9)), 
and then aim at obtaining finite limits of these values "per unit area" when the 
region increases to the whole plane. We shall not investigate the problem of the 
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thermodynamic limit in full generality and restrict attention to the following cases. 
The modified space JC± is a Hilbert space of functions on a rectangle with sides 
(La;, Lj,), the modified operator has an orthonormal basis of eigenvectors: 



hi_'ipki = m'fpki , '^ki = J {kf-xY + {le.yY , ex = tt/L^ , Ey = ■njLy , (2.3) 



and the scope of the values of (k, I) is cither 

(D) NxN, or (N) (N x N) U ({0} x N) U (N x {0}) . (2.4) 
For the scalar field the following choices lead to these two cases: 

(D) /C_L = L'^{{-Lx/2, Lx/2) X {-Ly/2, Ly/2),dxdy) and -h]_ is the two-dimen- 
sional Laplace operator with the Dirichlet boundary conditions; 

(N) IC± is the orthogonal complement of constant functions in the space 

L'^{{~Lx/2, Lx/2) X {—Ly/2, Ly/2),dx dy), and —h\ is the two-dimensional 
Laplace operator with the Neumann boundary conditions. 

In the case of Neumann boundary conditions we exclude constants, as they would 
lead to more singular case. Related, but slightly different choices, also lead to the 
cases (D) and (N) for the electromagnetic field (see Section [7|). For this field the 
absence of constants appears naturally. 

With these modifications we can now introduce external conditions. As a re- 
sult the z-motion Hamiltonian is replaced by a new positive operator hza, and 
in analogy to Eq. H2.1|l we have 

Viha)=V{h^(E)id)nV{id®h,a), ha = ^/{hi_ ® id)2 + (id®h^af ■ (2.5) 

A particular choice of hza will be proposed later, but first we want to identify 
all models fitting into the scheme, i.e. all perturbations which are admissible in 
the sense described in jlj. A preliminary condition for this is that the symplectic 
mapping La determined by h and ha as in Section 1-5 be bounded, with bounded 
inverse - condition (i) in that section. In the rest of the present section we find an 
equivalent form of this condition in terms of the z-motion Hamiltonians hz and 
hza, and express quantities Ma and £a in terms of these operators. In the next two 
sections the conditions for finiteness of A/'a and £a and their infinite plane limits 
per area are obtained. 

First of all we observe that each of the subspaces "^fci ® A^z is invariant under 
both operators h and ha- If we fix the basis {ipki} (we give up the freedom of phase 
factor multiplication) then each of the subspaces ipki ® ICz is naturally unitarily 
isomorphic to JCz- It is then easy to see that if we denote for u > 




V{h{u)) = V{hz) , 

V{ha{u))=V{hza) 



h(u) — \/h1 + uid , 



(2.6) 
(2.7) 



haiu) = \J hla + uid 
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then also 

Vih'^'iu))^V{hl/'), V(hl'\u))^V{h\i^) (2.8) 
and one can use the foUowing identifications 

/C = 0(/C.).,, /,V2^0;,i/2(,2^)^ h\/^^^hT{el)^ (2.9) 

where {K,z)ki are identical copies of /C^. The spectrum of both h^/'^{u) and ha^^{u) 
is contained in (u^/^,cx)), so both operators h^^^^ and ha^^"^ are bounded, with 
\\h^'^'\\,\\ha'^'\\ < (minfc,{efe,})-i/2. 

After these preliminaries it is now easy to show that the following conditions 
are equivalent: 

(i) The operators h and ha (cf. (|2.1|1 and (|2.5|1 resp.) satisfy the conditions 

= I?(/i±i/2) ^ (2.10) 
Ba = h^^^h^^^^ and B^^ extend to bounded operators in JC . (2.11) 

(i)' The operators and h^a satisfy the conditions 

Vihli^)^V{hl/^), (2.12) 
Ba{u) = h]J'^{u)h~^/'^{u) and Ba^{u) are uniformly 
bounded on each set u £ {v, oo) , v > Q . 

Note that the condition (i) is identical with the condition (ii) of Section 1-5, 
which is equivalent to the symplectic mapping La being bounded together with 
its inverse (condition (i) in Sec. T5). 

To prove the equivalence suppose first that (i)' holds. Then 

Ba = @Ba{eli), (2.14) 

M 

SO < max\\B^\eli)\\ < const. Then /i^^ = Bah^^^ and h^/^ = B-^h\'^ , 

ki 

so ft.^/^ and /i^^ have a common domain. As ft,^^/^ and ha^^^ are bounded, this 
ends the proof of (i). Conversely, suppose that (i) is true. Then the restrictions of 
the domains of /i^/^ and h]/^ to {ICz)ki must be equal, which implies Eq. (|2.12|) . 
Condition (|2.11() implies that Ba^{e1i) are uniformly bounded. But for u,w > 

Ba{u) = h]/^{u)h-^/\w)Ba{w)h^/\w)h~^'\u) , (2.15) 

so by spectral analysis for h and ha one finds 

115^(^)11 < max (^)'^'l \\B^\w)\\ . (2.16) 
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Setting here w = and u in the interval between e|; and a neighbouring value 
of e^, or between zero and e|;, if the latter is the minimal value, one shows that 
the condition (|2.13|) is satisfied. 

The equivalence shows that if (i) is satisfied for some L^^Ly, then it is true 
for all finite values of these parameters. 

If the equivalent conditions (i) and (i)' are satisfied then the operators Ta 
and Sa defining the Bogoliubov transformation are bounded, and given by 

"^a = Ta{eli) , Sa = Saiikl) 1 {2-17) 

kl kl 

where 

Taiu)^^[B~\u) + B:iu)], Sa{u)^^[B-\u)~B:iu)]K, (2.18) 

and K is the operator of complex conjugation. The quantities A/'a (particle number) 
and £a (Casimir energy), which have to be investigated, may be now written as 

Ma^TT[SaS:]^J2^T[Na{eli)] , 

kl (2.19) 
where Na{u) ^ Sa{u)S*{u) , 

£a^Tr[h'/'SaS:h'^']^Y.^r[Eaiel,)], 

kl (2.20) 
where Ea{u) = h^^^{u)Sa{u)S*{u)h^/^{u) , 

where we do not know yet whether these expressions are finite (but they are 
unambiguously defined as positive numbers, finite or not, as the operators are 
positive). 



3 The necessary and sufficient condition for the 
unitary equivalence of representations and the 
existence of finite Hmit Ua = hm Afa/LxLy 

LxiLy — >oo 

We now add the condition: 

(ii) The ground state representations determined by h and ha are unitarily equiv- 
alent, and the infinite plane limit of particle number per area is well de- 
fined, i.e. 



Ma < OO for all Lx, Ly , 3 lim — — ^ < oo . (3-f ) 

L^.Ly^ca LxLy 
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We show in this section that the conditions (i) and (ii) are satisfied if, and 
only if: 

(Eq) The operators and h^a have a common domain and hza — hz extends to 
a Hilbert-Schmidt operator on /Cz, i.e. 

Mh.a - K f < oo . (3.2) 



In this case 

Ha = hm — ^ = Tiihza - hzf . (3.3) 

In the case of Dirichlet boundary conditions for h± (case (D) in Eq. ^ZA\ ) the limit 
in Eqs. H3.1|l and (|3.3|) can be taken in arbitrary way, but in the case of Neumann 
conditions (case (N) in Eq. (|2.4|l ) one has to keep 1/M < L^/Ly < M for some 
arbitrary, but fixed M > 1 (which means that the limit takes place over the values 
of {Lx,Ly) within a conic proper subset of (0, cx)) x (0,oo)). 



3.1 (Eq) ^ (i) and (ii) 

Suppose first that the condition (Eq) is true, and denote c = \\hza — hz\\ - Then for 
£ ^{hza) = T^{hz) we have 

\\ha{u)^f - WhzaVf+uyf < iWhzifiW + cy\\f + U\\ipf 

< ('l + ^ + ^V||/i,^||2+^||^||2)= ('l + ^ + £!'j||/,(„)^||2. 

\ \/u u J \ u J 

Now we use the monotonicity of the square root: if A and B are positive operators, 
V{A) C V{B) and < PV'H for all ip G V{A), then also V{A^/^) C V{B^/^) 

and \\B^/^iP\\ < \\A^/^ip\\ for aU V' S V{A^/^) (see 0, proof of Thm.X.18). Thus 
we get 

This implies ||-Ba(w)|| < (l + cj ^fu + (? ju)^^'^ . Changing the roles of h[u) and 
ha{u) one exhausts the condition (i)'. 

We introduce a bounded selfadjoint operator on ICz- 

Aa{u) = ha{u) - h{u) (3.4) 

(boundedness by ||/iq(u) — /ia(0)|| < y/u, \\h{u) — h{0)\\ < ^/u and the boundedness 
of hza — hz). In the further course of the proof we shall need several theorems on 
Hilbert-Schmidt properties of operators like Aa(w) and Sa{u). As their demon- 
stration is somewhat more technical we shift their derivation to Appendix and 
in the main text use the results. We also remind the reader the following facts on 
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operators in Hilbert space: if B is bounded and H is Hilbert-Schmidt, then BH 
and HB are Hilbert-Schmidt, and if in addition K is also Hilbert-Schmidt, then 

Tt{BHK) = Ti{HKB) = Tt{KBH) . 

Furthermore, if B is in addition positive, then 

Tr {VbHH*Vb) = Ti{H*BH) < \\B\\ Tt{H*H) . 

We use these properties extensively in what follows. 

We have to show that the expression (|2.19|) satisfies condition (ii), and 
Eg. (|3.3() holds. We have assumed that Aa(0) is a HS operator. It follows then 
from the result (iv)B of AppendixlBlthat for all w > the operators Aa{u) are also 
Hilbert-Schmidt, and the function (0,c») 9 m i— *■ Tr A^(u) is continuous on its do- 
main and continuously differentiable on (0, oo), decreasing, and lim Tr A^(u) = 0. 

u — >oo 

Moreover, the differentiation on u may be carried out by formally pulling the 
derivative under the trace sign and differentiating formally the w-dependent oper- 
ators. Using the formal rule 

2^Aa{u) = -h-\u)Aa{u)h-\u) = - h-\u) Aa{u)h-\u) (3.5) 
au 

and the identity 

Sain) - -h-'^'{u)Aa{u)h-'/'{u)K (3.6) 
following from the formula H2.18|l one finds 

2^- TrAliu) = - Tr [h-\u)Aa{u)h-\u)AJu)] 

- Tr [Aaiu)h-\u)Aaiu)h-\u)] = -STvNaiu) . 

It is now clear that the function (0,oo) 3 u Tr Na{u) is continuous, positive, 
and 

/ TTNa{u)du= -Tt Alio) <oo. (3.8) 
Jo 4 

Again by the result (iv)B of Appendix IbI the function (0, oo) 3 u t-^ Tr Na{u) 
is decreasing. With the use of formulas (|A.14|1 and (|A.2I) of Appendix ^ this is 
sufficient to conclude that in the case (D) of Eg. (|2.4II the condition H3.1|) . and 
Eq. (|3.3(l are satisfied. In the case (N) of Eq. 1)2. 4|l J\fa is bigger by 

C30 

AA(^) - AA(^) = e^ey ^ Tr Na{{leyf) + {x ^ y} . 



1=1 



If one keeps 1/M < ex/^y < M then 



oo 



AAW _ A^i^) < Mel E iVa((/e,)^) + Mel E Tr Na{{lex f) 



1=1 1=1 
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which by (|3.8|) and the properties ljA.5|l and (|A.15|I of Appendix ^ is finite and 
tends to zero for ex,^y ^ 0. Thus in both cases (D) and (N) the condition (ii) and 
the equation H3.3|l follow from the criterion H3.2|l . 

3.2 (i) and (ii) (Eq) 

We now turn to the proof that the criterion (Eq) follows from (i) and (ii). This 
fact is of interest, as it shows that in the context adopted in this paper there is no 
escape from this condition. Suppose that the first of conditions (|3.1|) is satisfied. By 
the result (iv)B of Appendix IbI if Tr Na{u) is finite for any u > 0, then it is finite 
for all u e (0, +oo), and the function (0,+oo) 9 m i-^ Ti Na{u) is continuously 
differentiable and decreasing. Thus if Ma is finite for all L^^Ly., then by the first 
of inequalities in (|A.14p we have for each v > 0: 



poo 

/ Tr Na (u) du < oo . 

J V 



(3.9) 



Let now = P{o,p){hz) be the projection operator from the spectral fam- 
ily of the operator hz, projecting onto the spectral interval (0,(3), and simi- 
larly — P{Q.-y){hza) for the Operator hza- Then for each m > the expression 
Pf3Aa{u)Pa'y defines a bounded operator (although wc do not know yet whether 
Aa{u) alone makes sense), and 

Tr[PpAa{u)PajAa{u)Pp] 



< 4V(/32 + U)ij^ + U) TT[PpSaiu)Pa^S:iu)Pf3] 

<A^/(p^ + u){j^ + u)TrNa(u)^0 {u ^ oo) , (3.10) 



the last line by the property l|A.4p (note that here and in what follows all operators 
under the trace sign are positive). It follows from the result (iv)B of Appendix IBI 
that the function on the l.h. side may be differentiated as that in Eq. H3.7|) , so one 
finds 

^ TT[PpAa{u)Pa-,Aaiu)Pfj] - -4 Tr [P^j^a (u)P,^ 5: (u)P^] . (3.11) 

Integrating this identity on (i;,C5o), i; > 0, one has 

Tr[PpA4v)Pa^Aa{v)Pp] ^ 4 TT[PpSa{u)Pa^S:{u)Pp] du . (3.12) 

J V 

The integrand on the r.h. side is bounded by TvNa{u)^ so the double hmit 
/?,7 ^ oo exists. Thus 

rOO 

TiAl{v)=A TiNa{u)du. (3.13) 

J V 
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Now we let the second of the assumptions in (|3.1|) come into play. By formula 
ljA.2|l this implies that the r.h. side of the last formula has a finite limit for v ^ 0. 
As a result we have 

/>oo 

TtAI{v)<4: Tt Na{u) du < oo . (3.14) 
Jq 

By the result (iv)B this implies that Aa(0) is HS, which completes the proof that 
H3.2|l is a necessary condition. 

4 The necessary and sufficient condition for finite 
backreaction energy and its limit 

£a = lim Ea/LxLy 

LxjLy — yoc 

We now add the energy condition: 

(iii) The backreaction energy £a is finite for all L^, Ly, and the limit of Casimir 
energy per area is well defined for infinite plane limit: 

£a < oo , 3 lim < oo . (4.1) 

We show in this section that if the conditions (i) and (ii) (or criterion (Eq)) 
are satisfied, then (iii) is fulfilled if, and only if: 

1/2 

(En) The operator {hza — hz)hz is a Hilbert-Schniidt operator, that is 

Tr[{h,a - K)hz{hza - K)] < oo . (4.2) 
If this condition is satisfied, then 

£a= lim -^ = ^TT[{hza-K){2K+hza){Ka~K)]. (4.3) 
L^,Lj,— ►oo Lix-Ly Z47r 

The limit with respect to L^, Ly is specified as in (|3.3|l . 
4.1 (Eq) and (En) (iii) 

The proof is closely analogous to the one in the last section. Suppose that cri- 
terions (Eq) and (En) are satisfied. Wc have to show that the expression H2.20|l 
satisfies the conditions (|4.1|) and Eq. H4.3I) . We have assumed that Tr[Ao(0)^] < oo 
and Tr[Aa(0)/i(0)AQ(0)] < oo. It follows then from the result (v)b of AppendixEl 
that the functions u '—^ TT[Aa{u)h{u)Aa{u)] and u ^—^ TT[Aa{u)ha{u)Aa{u)] are 
continuous on (0, oo) and continuously differentiable on (0,oo), and tend to zero 
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for u — + oo. Also, the formal differentiatfon with respect to u as in Eqs. (|3.5|) . (|3.7|l . 
yields the correct result. A direct calculation then yields 

^ -^[2TT[Aa{u)hiu)Aa{u)] +Tr[Aa{u)ha{u)Aa{u)]} = -Tv Eaiu) , (4.4) 
and the integration leads to 



TvEaiu) du=- Tr[A,(0)(2/i(0) + /ia(0))A„( 
6 



< oo . 



(4.5) 



By the result (iv)B the function TvEa{u) is non-increasing. By a reasoning analo- 
gous to that following Eq. H3.8I) one shows that fa^"* — S'if^ vanishes in the limit. 
One concludes that the conditions (|4.1|) and the equation (|4.3f) are satisfied. 

4.2 (Eq) and (iii) ^ (En) 

Conversely, suppose that h^a — /i^ is a Hilbert-Schmidt operator and is finite. 
Similarly as in the case of T^r Na{u) it follows from the result (iv)B of Appendix lEl 
that if Tri?a(u) is finite for any u > 0, then it is finite for all u e (0, -l-oo), and the 
function (0, -|-oo) 9 u i-^ Ti Ea{u) is continuously differentiable and decreasing. 
Thus if £a is finite, then for each w > 



Tr Ea iu) du < oo . 



Therefore 



i:i[PpAa{u)Pa^ha{u)Pa-tAa{u)Pp] 

< (7^ + u) Tr \PpAa{u)^^Aa{u)Pp] < (7^ -|- u) Tr Ea{u) < 00 , 



(4.6) 



(4.7) 



Tv[Pa-Aa{u)PfMu)P(iAa{u)Pa 

< y/{f3^ +u){-f^ +u)Tt 



^^-^ Aa{u)PpAa{u)^'"' 



(4.8) 



< ^iP^+u){-f^ + u)TvEa{u) < 00, 



and by HA.4|I both functions tend to zero for u ^ 00. Differentiation is again 
allowed and gives 



^-^^2TT[Pa^Aaiu)Pi3h{u)Pi3Aaiu)Pa^]+TT[Pi3Aa{u)Pajhaiu)Pa^Aaiu)Pp]'j 

= - Tr 



F^A,(z.)-^A„(^)P^ . (4.9) 

ha{u) '^i 
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Now we take into the consideration the second condition in (|4.1|) and go over the 
steps analogous to those following Eq. I|3.11|l . This gives 



/•oo 



Tr[Aa{v)h{v)Aa{v)] < 3 / Ti Ea{u) du < oo . 



(4.10) 



Jo 



This implies, by (v)b, that Aa{0)h^^^{0) is also HS, which ends the proof. 



5 Modified Dirichlet and Neumann conditions 

From now on we can restrict attention to the z-motion dynamics, and specify 
a class of models. We set ICz = L^(M), — —d^, where the unique selfadjoint 
extension of the second derivative defined for S{M.) is meant. We want the external 
conditions to be some modification of the strict boundary conditions enforced on 
the planes z = ±6, where 2b = a. The Hamiltonian of the z-motion for such 
conditions is determined by (hf^)'^ = —{8^)3, where B stands for Dirichlet or 
Neumann conditions at z = ±6. Let F and G be real, non-negative, bounded, 
measurable functions on (0,oo). We postulate a class of models by setting 



We assume that for all m > we have F{u) < u and G{u) < 1, which guarantees 
the positivity of hza- The strict boundary conditions are recovered by formally 
setting F{u) = u and G{u) = 1 for all u > 0. Our intention is to keep these 
equalities for u "not to large" , and modify F and G so that for u 00 F remains 
bounded, and G (possibly) tends to zero. This seems to model correctly the idea 
that the boundaries should become transparent for very energetic particles. We 
show in this section the following: 

(Mod) Let the functions F E C^((0, 00)) and G E C((0,oo)) satisfy the estimates 

< F{p) < const , F{p) <p, < G{p) < 1 , 

G(p) < const(p + 1)-" , |F(2)(p)| < const(p+ l)-2-7 (5.2) 

for some a, 7 G (0, 1) and such that 4q; + 27 > 1 . 

Then hza — hz, where hza is the modified z- motion operator given by (|5.1|) . 
satisfies conditions (Eq) and (En) . In this case hza — hz is an integral operator 



iij, [hza - hzW) = J i^{p) <''=(p,p') V^'(p') dpdp' , (5.3) 

with K^'^{p,p') e L2(]R2) and ^/\f\K^-^ ip,p') e L^{M.^). The explicit form 



hza = hz + G(hz) [F(/ifj - F(hz)\ G{hz) . 



(5.1) 
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of K^'"^ is obtained in (|5. 231 - 15. 2511 below. The Casimir energy is given by 

= i{ / dpdp' (5.4) 
+ j K^'^ip,p')K^'^ip',p")K^-^ip",p)dpdp'dp"Y 

Moreover, let for /i > 0: 

F^{p)^^iFip/^l), G^{p) = G{p/ii), (5.5) 

and denote by the Casimir energy for the model with G replaced with 
F^,G,,. Then 

e'^a = ^J?e^a ■ (5.6) 

As explained in the opening paragraph of this section, we believe that (|5.1() , with F 
and G in suitable classes, defines a class of models reasonably imitating boundary 
conditions for low energies, while giving the means for "softening" of the walls. 
More specific of the assumptions on F and G are technical. (Mod) then shows 
that indeed those models are not only admissible for the discussion of adiabatic 
backreaction as explained in earlier sections, but also susceptible to a rigorous 
treatment. Low energy behaviour F{p) ~ p and G{p) ^ 1 for "small" p, important 
for the above interpretation, is not needed in (Mod). Once these assumptions 
are added in the next section, the scaling properties described in 1)5. 5|l and 1)5. 6|l 
provide means for investigation of the sharp boundary limit. 



5.1 Spectral representation of h 



B 
za 



The spectral representation for hz is obtained by Fourier transforming i— > -0 
{ip,i^ G L^(R)), so that in the new representation hz acts by hzipip) — |p|V'(p)j 
where the conventions for the transform are defined by 

ipip) = / tPiz)e-'P''dz . (5.7) 



/27r 

To prove (Mod) it is sufficient to show that in this representation F{hf^) — F[hz) 
is an integral operator with a kernel K^{p,p') such that both functions 
Gi\p\)K^{p,p')G{\p'\) and y^\G{\p\)K^ {p,p')Gi\p'\) are in L^{R^). 

We start by identifying the spectral representation for hf^. It will be conve- 
nient to denote 

I — 1 Dirichlet case , 
= < (5.8) 
+ 1 Neumann case , 
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and introduce 

Ll = {xeL\R)\x{-z) = Sx{z)}, lj = {x&l\^)\x{-k) = Sx{k)}. (5.9) 

For ijj e i^(M) (in the initial position representation) we denote 

i;.{z)=e{-z-b)ilj{z), Mz) = 0{z + b)e{-z + b)i,{z) , i;+{z) = e{z-b)iP{z). 

(5.10) 

Then the following rule defines a unitary transformation: 

L\R) 9 V ^ * = (V^ ,V^,V?) e e /| ® 

— nr / ^ ^ \ TT (5-11) 

^oik) = ^ ^(i'Mke) + Si-'^M-ke)) , e=-, 

V?(r) = i= [e+"'-ir+{r)+Se-">^ir+i-r)) , 

where , V^o , V'^ ) arc Fourier transforms of (■(/>_, -i/^Oj ''/'+) respectively. In this 
new representation operator hf^ acts on ip±{r) as multiplication by \r\, and on 

tpf,{k) as multiplication by \k\e. This is easily seen by using an overcomplete set of 
eigenfunctions in the three regions 



^±ir,z) = -J_(e-(^T6) (^Tb)) ^6ip±{-r,z), z ^ ±b 

Zy/TT 

Mk,z) = -^{i-'^e"''' +Si''e-"''') = dM-k,z) , z G (-6,6) 
a 



(5.12) 



normalized by 

j (p±{r, z)ip±{r' , z)dz = ^[6{r - r') + S 5{r + r')] , 
J ipo{k,z)ipo{k',z)dz = ^[Sk,k' +^^fc,-fe']- 

We also note for later use that for tp G S{M.) by standard Fourier methods 

— If smb{p-r) - 

tpo{r) = - / tpip) dp , 

IT J p — r 

2n J p — r — lO 2-K J p — r — lO 

ZTT J p — r + iO ZTTjp — r + iO 



(5.13) 



(5.14) 
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5.2 The integral kernel of hza — hz 

Let V'jV'' G 5(M). Using the spectral representation (|5.11l) one finds 



(v,F(/ifjv')= pi\r\) 



dr 



+ S / F{\r\) 



'■V'-(-r-)V'- W + e""'V'+(-f)V'+W 



(5.15) 



For further evaluation we shall need results of the Appendix [CI We also use the 
notation introduced at the beginning of that appendix. 

We represent the terms in brackets under the integral sign in the first line by 



tjj-{r)ip'_{r) + tp+{r)ip'^{r) ='ip{r)'ip'{r) — 4>o{r)%l;' {r) 



(5.16) 



and use representation given in H5.14|l . Then for the integration of the terms in 
the second line of the last identity we can use identity 1C.21|I . We get 



F{\r\) tlj^{r)tp'_{r) +ip+(r)'il;'{r) dr = F(\r\)ilj{r)il;'{r) dr 



+ 



1 

2^ 



\Kp-p')]xf{p,p') 



(5.17) 



o / ML- vi^i^ ■ - i/j ip{p)^''{p)dpdp' 
For the integral in the second line in H5.15|) we again use l|C.21|) to obtain 



F{\r\) e-'"'V-(-r)V'^(r) +e"'^V+(-?-)V'VW 
1 



dr 



2ti 



2 cos[h{p ~ p')]xf{p,~p') 



(5.18) 



'4'{p)iIj'(p') dpdp' . 



2tt{p + p') 

To evaluate the sums in the last line of Eq. H5.15|l we represent V'o as in 1)5.1411 and 
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use identities (|C.24|I and (|C.25|) . which results in 

fcez 



cos[b{p- p')]xfAp,p') 



cos[b{p + p')][xfAp,p') - XfMp^p')] 



(5.19) 



ipipW {p') dpdp' , 



6^i^(|fc|e)(-lfM-fce)V^(fce) 



^ cos[6(p - p')]xfAp^ -P') 



^ cos[&(p + p')][XKe(p, -p') - XfMp^ ~p')] 



(5.20) 



i^{p)4''ip') dpdp' . 



Finally, we note that 



F{K)i,') = j i^(|r|)V'(r)^'(r) dr . (5.21) 
Setting now expressions IjS. 171 - 1^?^ into H5.15|l and subtracting H5.21|l we find 

(V, [^^(/ifj - F{h,)W) = j J{p)K^{p,p') ^'{p') dpdp' , (5.22) 
and as a consequence obtain Eq. H5.3|l . where the kernels are given by 



K^{p,p') = cos[6(p -p')][fcF + kFA(p,p') 

+ COs[b{p + p')] [kF.e - kF,2e] {P, p') 

= cos[h{p-p')\[kF + k*pj{p,p') + cos[6(p + p')][fcF.e - k*p22{p,p') 

+ 2sm{bp) sm{bp')k°p^{p,p') , 



with 

kF{p,p') = Xf{p,p') - Sxf{p,-p') 

kF,e{P'P') = ^ XfAp^p') - ^XfAP' -P') 
and similarly for fc^^, k*p and 

<^«(p,p') = G(b|X(p,p')G(b'|) 



(5.23) 



(5.24) 



(5.25) 
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(all denotations as at the beginning of Appendix [nj. Using the second formula in 
H5.23|l and looking at the estimates (|C.16I - IC.19I) one finds that both functions 
K^''^{p,p') and y/\p\K^''^{p,p') are square-integrable. In more detail: square in- 
tegrability on the square {p,p') E (—1,1)^ follows rather immediately from the 
estimates (IC.16l - rC.18|l (the term containing kp^ has to be analysed jointly with 
the sine functions - see H5.23|l \ For the square integrability on \ (—1, 1)^ one 
uses HC.19|) and H5.2|) ; it is obvious that in this case the integration of the estimates 
can be harmlessly extended to R^, and then by symmetry narrowed to (0,oo)^. 
In this way one finds 



[ [i+H]K'<=(p,p')rrfprfp' 

Jm-^\{-i,i)^ 



where for the integration over p > p' > the variables have been swapped, and the 
condition 4a + 27 > 1 has been taken into account (cf. (|5.2|) '). The formula H5.4|l 
for the Casimir energy is now easily obtained (cf. Eq. H4.3|l 'l. Finally, to obtain the 
scaling behaviour it is easily checked that Ka'^' '^{p,p') = K^^{p/ ^,p' / which 
upon substitution into H5.4|l gives the desired result. 

We end this section with remarks on the relation of the present work to 
that of pp. The models considered in the cited work are obtained by choosing 
the Dirichlet boundary conditions for h^^ and setting G = 1. One may thus ask 
whether it is possible to set G = \ also in the case of Neumann conditions, which 
would alow to bring the analysis to simpler terms of It turns out that it 
is not, and this posed the major difficulty in extending the program to these 
boundary conditions. The estimates (|C.16I - rC.19|l guarantee that K^{p,p') is 
square-integrable, but are insufficient for the square-integrability of ^/\p\K^ {p,p') 
(this is easily seen by putting a = in H5.26|l '). However, if one takes into account 
that the functions x enter the kernel only through the combinations H5.24|l . then it 
turns out that the terms in x's decaying most slowly for \p\, \p'\ ~* 00 subtract in 
the case of Dirichlet, but add in the case of Neumann conditions, and the hmction 
^y\p\K^ {p,p') is square-integrable in the former, but not in the latter case. 



6 Asymptotic expansion of the Casimir energy 

In this section we show how to expand the Casimir energy of the class of models 
defined in the previous section into inverse powers of the distance of plates a. 
We shall assume that F and G are as for strict boundary conditions in some 
neighbourhood of zero, F{p) = p and G{p) — 1 for small p. This assumption 
is somewhat stronger then the most economic one, which demands only that a 
few derivatives of those functions in p = agree with those of the above special 
functions. We choose the simpler version to shorten the proofs. Also, we assume 
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the differentiability to arbitrary order, wfiicfi spares us tedious track- keeping. We 
note that to take /i arbitrarily large in H5.5|l and (|5.6|l physically means to approach 
strict boundary conditions. The asymptotic expansion combined with the scaling 
determines the behaviour of the Casimir energy in that limit. 
Below we prove the following. 

(Asym) Let the functions F,G E C°°[{0, oo)) satisfy the conditions 

< F{p) < const , F{p) <p, < G{p) < 1 , 

<const(n)(p+l)-"-^ n= 1,2,... , 
\G'^^\p)\ < const(fc) (p + l)-*^-" , fc = 0,1,2,... 
for some a, 7 G (0, 1) and such that 2(a + 7) > 1 , 



(6.1) 



and let moreover F{p) — p and G{p) = 1 in some neighbourhood of zero. 
Then 

2 

ea=e^ + ( - + ^4,a , (6.2) 

where Eoo and c are model {F and G)-dependent constants and 



. . const 
\£aA < — 5- 



da 



const 

<-:^. (6.3) 



The limit value £00 is twice the Casimir energy eo for the configuration of 
one single plate, 

£00 = 2£o . (6.4) 

The scaling limit /i ^ 00 yields 

lim (£(^ - £oom' - (i±^)4/^) = -t£^, ■ (6-5) 

Constants £00 and c are given in H6.9|l . (|6.32|) . H6.52|l and 1)6. 53|) below. 

As already mentioned above we do not strive to optimize assumptions on functions 
F and G, which could be substantially weakened, but rather want to simplify the 
proof. The assumptions on derivatives essentially mean that the functions do not 
have oscillatory behaviour; the stronger (then in (Mod)) bound on parameters a 
and 7 will be needed in our expansion procedure. 

The usual interpretation of formula 1)6. 5|l in the Dirichlct case would be that 
apart from an "unimportant constant" the limit of strict boundary conditions 
gives the original Casimir expression. However, this is not true for the Neumann 
conditions. Moreover, the constant £00 does have physical meaning: this is the en- 
ergy needed for the creation of the configuration of the field surrounding infinitely 
separated plates. Relation 1)6.4)1 shows that configurations surrounding each of the 
plates become independent in the large separation limit (at least as far as energy 
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is concerned). We also note, although we do not present explicit calculations, that 
in the sharp boundary limit also the particle number per area diverges. 

One could ask how generic are the lessons drawn from (Asym). We think 
that the result whose generality is least prone to doubt is the nonexistence of sharp 
boundaries limit. We support this by two observations. First, we recall that the free 
quantum field and the field in presence of sharp boundaries are even not described 
by the same algebra of observables, not to mention their representations (this 
was pointed out in P). Second, we learn that quantities which are responsible for 
comparability of representations and for backreaction, and Sa, become infinite. 
Now, note that these are positive quantities, having physical interpretation. These 
facts indicate very strongly to our conclusion. Similar views on the unphysical 
nature of sharp boundary conditions may be found in literature, but we are not 
aware of a rigorous discussion based on the investigation of the algebraic structure 
of the theory and modification of dynamics (see also discussion and bibliography 
in D). 

Also, we note that our calculation offers a rigorous derivation of the original 
Casimir term — 7r^/1440a^ in the backreaction energy and confirms its universality. 
On the other hand the derivation shows, that it cannot stand alone: the energy is 
the expectation value in certain state of a positive operator, so it must be positive. 
Even if one adds a positive, constant (a-independent) term to the Casimir term 
one does not get a positive expression for each separation value a. 

The appearance of a term quadratic in 1/a in the Neumann case stands 
probably more open to debate. However, we want to point out that our rigorous 
analysis should rather be viewed as giving rise to a converse problem: are there 
any models for the dynamics hza, for which the quadratic term would be absent?; 
what would distinguish them? Within our class of models there is no general reason 
for c to vanish, and at this point we are not aware whether the condition c ~ 
has any solutions in that class. In fact, a partial result to the opposite in certain 
circumstances may be shown (see the end of this section), for brevity we simplify 
assumptions on functions F and G to a more special class then needed. Suppose 
that in addition to the assumptions of (Asym) we have G{p) = f{p — r) for p > 0, 
where r is a positive parameter and f{x) is a smoothed step function equal to 1 for 
X < —K and to for x > k, with some k > 0. Then for sufficiently large r constant 
c is positive. Note that increasing r moves us towards better approximation of 
sharp boundaries. Further discussion of the physical meaning of our results will be 
found at the end of Section 

We shall assume in the proof of (Asym) that F{p) — p and G{p) = 1 for 
p G (0, 1). This does not restrict generality: each case of functions satisfying the 
assumptions of (Asym) may be brought to this more restrictive case by rescal- 
ing H5.5|l . 
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6.1 Asymptotic value of the Casimir energy 

We first consider tlie Casimir energy for the z-motion dynamics of modified single 
plate. Let /ifo be determined by (/ifo)^ = ~(^^/ (^-^)^) b with Dirichlet or Neumann 
conditions on the plane z = 0, and define h^o as in Eq. (|5.1|) . The discussion of 
the last section simplifies greatly in that case and shows that hzo — is again an 
integral operator in the momentum representation, with the kernel 

- G{\p\)kF{p,p')G{\p'\) , (6.6) 

with kp as defined by (|5.24l) . The Casimir energy per unit surface eo for this 
operator (formula (|5.4I) with this new kernel) is the energy one needs to create the 
configuration of the field in the lowest stationary state of the modified Hamiltonian, 
as discussed in [J. 

We now consider the limit value £oo in our model of modified par- 
allel plates. For a oo, i.e. e 0, one has for p 7^ the point-wise 
limit limAi?_e(p) = Af{p)- This is easily shown with the use of for- 

mula (|A.12|I and the estimates (|C.12|I and (|C.13p . Then for p,p' ^ 0, 
p ^ p', there is \imx*F,eiP>P) = Xf{p,p), and MuixfAp^p) = 0. At the 

£—►0 ' e^O 

same time referring back to the estimations discussed towards the end 
of Section |S1 one realizes that the functions G'(|p|)[xj-c — Xf]{p,p')G{\p'\), 
^/\p\Gi\p\)[x*F,e - Xf]{p,p')G{\p'\), sin{bp)smibp')Gi\p\)x%/p,p')G{\p'\) and 
^y\p\ sm{bp) sm{bp')G{\p\)x% ^{p, p')G{\p'\) remain bounded in modulus by square- 
integrable functions independent of e (when it is small). Hence for e — > we 
have 

\\Gi\p\)[x*F,,-XFKp,p')Gi\p'\)\\L2^0, 

WVWlGmixh - XF](p,p')G(b'|)|L. - , 
||sin(6p) sin(6p')G(|p|)xL(P,P')G(b'|)||L2 ^0, 

||/Hsin(6p)sin(6p')G(b|)x^,.(P,p')G(b'|)||^. ^0, 

so for the purpose of calculating the limit a ^ 00 one can replace the kernel 
Ka'^(p,p') with cos[b{p - p')]2K^'^(p,p') in the formula (El- Using now the 
identities: 

cos2[6(p-p')] = \{l + cos[2b(p-p')]] , (6.7) 

cos[6(p — p')] cos[6(p' — p")] cos[6(p" — p)] 

1 ( (6.8) 
= - 1 f + cos[2fe(p - p')] + cos[25(p' - p")] + cos[26(p" - p)] | 
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together with the Riemann-Lebesgue lemma we find 
£00 = J 3\p\[K!>^''iP,P')?dpdp' 

+ f K^ ''{p,p')K;^ ''{p',p")K^'''{p",p) dpdp' dp" ^ = 2eo . (6.9) 



6.2 Decomposition = + ^^"^ 



We now turn to the derivation of further leading terms of the asymptotic expansion 
of Sa for a — > 00. It wiU be convenient to denote 

= ^(/^fa) - n/l.) (6.10) 

and 

g{u) = l-G\u). (6.11) 

We note that the operators and ^/hlG{hz)D^ are HS - this is shown by 
methods used towards the end of Section [3 for the estimation of kernels (and here 
for the second one of these operators we need the strengthened condition on a 
and 7, see (|6.1|l ). Using this notation and manipulating the operators under the 
trace sign we obtain 

Tr(/.,, - Kf = Tr )'(1 - igih.))] 

„ (6.12) 
+ 3Tr [{D^,fg{K)D^,g{K)] - Tr [D^ g{K)D^, g[K)D^, g{K)] , 

Tr [{Ka - K)K{Ka - K)] = Tr [V^G(/i,)(i?f )'g(/i,)0^] 

-TT[D^,G\K)KD^,g{K)\. 

We set these expressions into the formula for the Casimir energy, in addition we 
add and subtract the term Tr \D^F{hz)D^{l — 3.g(/iz))] , and write the result in 
the form 

= (6.14) 

+ sVKMhz ) {DlfG{K)^, ~ {Dl) V(/i, ) } , 

Tr [iD^ [F{K) - G^{K)h,] D^^giK) 
247^ ^ (6.16) 



/J 



3D 



^)'g{h^)D!g{h^) ~ D^g{K)D^g{h,)D!g{h,)} 



This decomposition of Sa is purely technical: different techniques will be now em- 
ployed for the calculation of the asymptotic expansion in each of the two cases ei'' 
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and Ea ■ However, it could help to observe that if it was possible to put G = I 
(as can, in fact, be done in the Dirichlet case - see remarks towards the end of 
Sectional), then the second term in H6.14|l would vanish (as then g = 0). More 
comments on this point will be found below. 

6.3 Integral representation of s^'^ 

By elementary manipulations one finds the identities 

(D^)' = - F{K)D^ - D^FiK) , (6.17) 

(D^f + D^F{h,)D^ 

3 2 2 (6-18) 

= D^ -Fih,)D^ -D^ F{h,)+F{K)D^F{K). 

Operators and D^^^ are HS, with the kernels in momentum representation 
G L'^{M?). The above two identities, when written in terms of kernels, 
take the form 

I K^{p,q)K^iq,p')dq = Kfip,p')-[Fi\p\)+F{\p'\)]K^ip,p'), (6.19) 



K[{p,q)K[{q,q')K^{q',p')dqdq'+ K[ {p, q)F{\q\)K^ (q,p') dq 
J J (6.20) 

= Krip,p') - [F{\p\) + F(b'l)] iff (p,p') + Fi\p\)K[{p,p')Fi\p'\) . 

These identities hold in the L^-sense. However, as both sides are continuous func- 
tions of {p,p') for p,p' ^ 0, they are also valid there point-wise, in particular also 
for p ~ p'. Using this fact for the calculation of Ij6.15|l one finds 
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[l-3.g(b|)]iff (p,p) 

-3[G'i\p\)\p\~Fi\p\) + 2Fi\p\)gi\p\)]Ktip,p) (6.21) 
+ 3G\\p\)Fi\p\)[Fi\p\)~2\p\]K^ip,p)\dp. 



We note that this expression engages only the diagonal values Ka{p,p) of the 
kernels. In fact, as part ea"* is concerned, one could bypass the calculation of the 
kernels and use a simpler technique employed in Ref. 1 . But once we have the 
kernels (needed for the calculation of ei"-*), it is convenient to use them also here. 
Looking at eq. (|5.23l) we see that 

Ka{p,p) = [kp + kF,e]{p,p) + cos{ap)[kF,e - kF,2e]{p,p) ■ (6.22) 
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Using the definitions HC.4I) and (|C.7p wc find xf (p, p) ~ (p) and xf {p, —p) = 
Ap{p)/p and then recalling (|5.24(l we have 



kF{p,p) 



1 

4^ 



1 



p-q p+q 



■[^^(kl)-^^(bl)]- 



jdg. (6.23) 



In a similar way one finds 



ke 



[F{\ke\)-F{\p\)]- 



2|p|fW(H) 
{q-p){q + p) 

2\p\FW{\p\) 



{he — p){ke +p) J 
(6.24) 

Also, analogous formulas for or replacing F in H6. 221 - 1^.2411 are valid. Setting 
all these formulas into Eq. H6.21|l one finds that the terms proportional to i^*^^-'(|p|) 
cancel out and one is left with 



.0) 



1 



247r3 
1 



p{q,p) dqdp 



247r3 
where 

p{q,p) = A|(g,p) 



/ I [1 + cos(ap)] e ^ p(fce, p) - cos(ap) 2e ^ p(2fce, p) | dp 



(6.25) 



^^p^ + ^q' 



F{\q\)[l-39{\p\)]~Fi\p\)+3G\\p\)\p\ 



(6.26) 



All integrals and sums are absolutely convergent. Both here and in further ex- 
pressions below involving p and its derivatives one uses the estimates found in 
Appendix mi to prove convergence. We do not go into easy but tedious details. 

6.4 Expansion of e^^ 

As p{q,p) is even in each of the variables, one can reduce integration and summa- 
tion to nonnegative values. We evaluate the sum e p{ke,p) for fixed p > with 

fcez 

the use of formula (|A.12p for n — 4: 

e^p(fce,p) = 2/ p{q,p)dq-e^B2P^^'°H0,p)-2e*Rp^i^,{p), (5.27) 



where 



1 1'°° 

i?P,4,.(p)-^y^ P^^^°Hq,P)b4{q/e)dq 



(6.28) 



and p''^''^^{0,p) is the limit value of p'^^'^'^ {q,p) for q \ 0. To obtain the other sum 
in the integrand of Eq. H6.25|l one needs only to replace e by 2e. The integrability 
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of |i?p^4_e(p)| at infinity is easily shown with the use of estimates to be found in 
Appendix C. On the other hand, these estimates do not guarantee the integrabihty 
for small p. However, this property will be obvious from what follows below. Setting 
the above expansion into Eq. (|6.25|) we have 



1 pOQ 2 pOQ 

4'^=Tt/ p{q,P)dqdp~-^ p(i.o)(0,p)[l-3cos(ap)]dp 
3- Jo 72. J, 

I - [1 + cos(ap)] e^RpA,e{p) + cos(ap) {2e)'^RpA,2eip)j dp . 



6tt^ 

The first term on the r.h. side is a-independent and it contributes to £oo- To eval- 
uate the second term we find explicitly 

p(i'")(0,p)^-3(lM)^(p), ^(p)^G2(p)^M2P_ZM>o. (6.30) 

Using this one finds 

3 /i'°)(0,p)[l-3cos(ap)]dp 
^^'^ (6.31) 



= (^)i2wX '^'^''^^^l .«(,)sin(Mdp|, 

where integration by parts in the second term has been performed. From H6.3U|I 
we find that ri{p) = 1 for p e (0, 1), so "q^^^p) = in this interval. It is easily 
seen by standard Fourier transform properties that the second integral in the last 
formula and its derivative with respect to a vanish faster than any inverse power 
of a for a — > oo (by repeated integration by parts and use of estimates (|6.1|l '). The 
first integral contributes only to the constant c in (|6.2|l . and yields 
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Next, to evaluate the second line on the r.h. side of Eq. H6.29|l we split 

RpA,€ = RprA,c + Rp^A.e, where 

p^ p^ 

p = p, + Ps, ps{q,p) = -, — + (6-33) 

[p+qY p + q 

For q,p e (0,1) we have p{q,p) — ps{q,p) + (1 — S)q/2, so pr is regular in a 
neighbourhood of zero. Thus using H6.28|l we find the estimate 



^ J I - [1 + cos(ap)] e''i?p^,4,e(p) + cos(ap) (2e)''i?p^,4,2e(p)| dp 



^ const, |^(4.o)(^^p)|^^^p<^. 



(6.34) 
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For the integrability of one uses the assumption (I6.1|) . relation H6.33|l 

and the results of Appendix C. Furthermore, using the relations 



f°° dcosjap) I dcosjap) 
J Rp,AAP)d-P = - J pRprAAP)dp 

1 

= / cos(ap)[pi?2^4 ^{p) + Rp^aAp)] dp , 

"Jo 



(6.35) 



1 

<!4,.(P)-^/ pi^''Hq,p)b,{q/e)dq, (6.36) 



= 5e3i?^„4..(p) + ^ qpi^'°Hq,p)b4{q/e)dq, 



(6.37) 



one also finds 



la 6^ / { ^ ^ cos(ap)] e^Rp^AAp) + cos(ap) (2e)''i?p,^4^2e(p)} 



const 



(6.38) 



dq dp < oo . 



Integrability of p\pr'^'^\q,p)\ and of q|pr^'''''(<;,p)| is shown as for \pr'^''^\q,p) \ above. 

The second term resulting from the split of p, that involving Rp^A,<^' '^^^ t)e 
explicitly evaluated. We have 



f^^Rp^A.Ap) = e^-Rp3,4,i(p/e) , 



(6.39) 



Rp.AAt) = i-t'dt + St')- (^dt-^)bAu)du. (6.40) 



With the use of formula (|A.8p with n = 4 and e = 1, applied to the function 
f{u) = l/{t + u) (for fixed t) we find 



= lim 

Af-+oo 



N , N 

du 



t + u 



k=0 

1 B2 _ W2{t) 

2t 2^2 - i2 



1 ^2 

i + fc 1 ' 2t 2*2 



(6.41) 



where we used formula IjD.lj) and notation introduced in IID.6p . and then 

Rp^AAt) ^ -tw2\t) + ("2 + ^)Mt) ■ (6.42) 
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Setting this into H6.39|l and using the identities (|D.8|I and (jD.9|) we find 
^ J { ^ [1 + cos(ap)] e'^Rp^^iAp) + cos(ap) (2e)'' i?p^, 4, 2e(p)| dp 

= TT^ / [1 + cos(7ri) - 8 cos(27rt)] \tw^^'> (t) - (2 + S)w(t)] dt (6.43) 
60^^ ./ 



1 



rC(4) = - 



167r2^' ' 1440 a3 ' 

To arrive at the last hne in this formula one has to take into account on the r.h. 
side of SdM that 



00 ^ 00 ^ 

6.5 Integral representation and expansion of e^"^ 

We now turn to ei' - the second term in (|6.14() . We use identity H6.17|l in formula 
H6.16|l and calculate the trace in the momentum representation, which yields 

4"^ = ^/ K[ip,p')K[\p',p)gi\p\)gi\p'\)dpdp' 

- ^ / [K[ip,p')]%{\p\){G'{\p'\)\p'\ - Fi\p'\) + 2F(b'|)5(b'|)} dpdp' 

K^{p,p')K^{p',p")K^{p",p)g{\p\)g{\p'\)g{\p"\)dpdp'dp". 



1 

24^ 



(6.44) 



This expression looks much more complicated than the integral formula H6.21(l for 
Ea^ (multiple integrals and full kernels Ka{p,p')). However, a substantial simplifi- 
cation of the asymptotic expansion will be due to the fact, that the integrations 
here stay away from the singularities of the kernels Ka{p,p'), which occur in the 
neighbourhood of p = or p' = (as for u £ (0, 1) there is G{u) = 1, g{u) = 
and F{u) = u according to our assumptions) and the integrands are, consequently, 
infinitely differentiable. Namely, products of kernels Ka contain products of cosine 
functions, which may be expressed as linear combinations of other cosine func- 
tions and, possibly, unity (as in Eas. (|6.7|l . 1)6.81) ). Those terms in the integrals 
above which contain one of the resulting cosines yield functions of a which van- 
ish faster than any inverse power of a for a — > 00. This is shown by standard 
Fourier methods: integrate by parts and use decay properties of the integrand 
and its derivatives - estimates (|6.1f) and HC.20|) : in the process the decay rate 
of the integrand increases. One also shows that the a-derivative of those terms 
has similar decay properties. Indeed, the derivative when acting on trigonometric 
functions does not change the decay; if it acts on a sum of the form /(fee) it 
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yields —{^/a)e'^[f{ke) + kef^^\ke)]. These facts are used for the proof of the last 
statement, but we do not go into further straightforward details. 

Thus it is sufficient to consider only those terms which take a constant from 
the product of cosines. One finds easily that the only products of cosines in the 
above integrals which do give constants are those in which cosines with "+" sign 
between the variables do not occur, or occur twice. However, we shall see below 
that these parts of the kernels K which multiply the cosines with the "+" sign are 
of order at least, so the quadratic terms of that type do not contribute to the 
orders < 3 in e in the energy, and to the orders < 4 in the force. Thus we shall be 
left only with the products H6.7|l and (|6.8|l . 

To expand the kernels Ka we use the definition (|C.5I) and the formula ljA.12|l 
to get 

^2 



bp 



where 



1 f°° 

Af'aAp) = Xf^\q,p)h^{q/e)dq, 
Setting these expressions into IjC.Sp and 15.24|l one finds 



1,2. 



kF,e(p,p') = kF{p,p') 



1 + 5 
127r^ pp' 



e^k 



FA. 



■i{p,p') , 



kF^AP^P') = kF2{p,p') 
where for i — 1,2: 

' AF^A,eij>)- AFiA,eijp') 



^F'A. 



.{p,p') 



1 

2^ 



e ^f^aAp^p') ' 

Af^aAp) 



■A 



F\A, 



dp') 



(6.45) 
(6.46) 

(6.47) 
(6.48) 
(6.49) 



p — p p + p 

Note that \p\, \p'\ > 1 for our present purposes, so there are no singularities. Then 
with the use of bounds similar to those in (|C.20|) one shows that the terms con- 
taining kFA,e indeed do not contribute in the given orders. Tlius omitting the last 
term in both formulas H6.48|l and taking into account the result of the discussion 
following Eq. (|6.44|) we find that disregarding terms higher than (l/a)"^ in ei"'' we 
can make the following replacements in H6.44|l (recall the formula ()5.23f) 'l: 

\2kF{p,p') 



Ki{p,p')Ki [p'p) 



[<{pyf 



2kF{p,p')- 



1 



127r^ pp' 

^2 n2 



2kp2{p\p) , 



12n-^ pp' 



K^{p,p')K^{p'y)K^{p",p) 



2kF{p,p')- 



1 + 5 
12tt^ pp' 



2kFip',p") 



1 + 5 
127r2 p'p" 



2kF{p",p)- 



1 

r27r2 p"p 



(6.50) 
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In this way we obtain 



c(") 

r^") = (iii) — + terms independent of a + Oia^^) , (6.51) 



where 



(ii) ^ f 1 ( n9{p)g{p ) ^ , I 

- =--r7r kp2{p,p) ; dpdp 

Utt Jq pp' 



kF{p,p) 



on Jq pi p' p' 

OTT Jo P P 



dpdp' (6.52) 



The constant c in 16.211 is now 

c = c(') + c('') . (6.53) 

This ends the proof of Eqs. H6.2|l and (|6.3|l . Equation (|6.5() is their simple conse- 
quence. 

FinaUy, we want to prove the statement on positivity of c in circumstances 
described towards the end of the discussion foUowing the formulation of (Asym) . 
We note that the term c*^'-* obtained in (|6.32|) is then a positive, increasing function 
of r tending to the limit value 

lim c(0(O^-i- r^(^)[^^;^(^)],,. (6.54) 

r^oo 247r Jo P 

At the same time one shows with the use of the present assumptions on G{p) and 
estimates on functions k that |c^")(r)| < const /r''', so 

lim c(")(r) = 0, (6.55) 

r — *oo 

which ends the proof. 



7 Electromagnetic field 

When dealing with a scalar field we have first formulated the free field model in 
terms of the value of the field, and its time-derivative, on a Cauchy hyperplane. 
Then we formulated the dynamics in presence of the external influence in terms of 
the same variables, and took over for the unperturbed energy of the field itself (at 
a given time) the old expression built with the use of them. In this way a simple 
model could be obtained in which the "sources" of the field were supplied by the 
field variables themselves. A similar procedure for the electromagnetic field needs 
some caution. The main problem lies in the fact that the dynamics of the field 
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is constrained, and one of the constraints depends on sources (Gauss law). The 
evolution is governed by a first order equation and one could think of the electric 
and magnetic parts of the field as the analogues of the scalar field variables at a 
given time, but then the initial data for the electric part are differently constrained 
in free and interacting case. Therefore one first has to solve constraints and identify 
independent field variables. The modification of the dynamics of these variables 
must then produce an imitation of the perfect conductor boundary conditions 
of the electromagnetic field. To achieve this, one has to adjust the choice of the 
independent fields to the geometry of the problem. We shall first describe our 
choice in the case of the classical free field restricted to the region x K, where 
is a compact convex region in the x-y plane, and M is the z-axis. The following 
informal discussion serves to motivate the precise formulation of the model, which 
starts after Eq. ^A^. 

For a vector A we denote by A± its part perpendicular to the z-axis, and by 
A its 2;-component. By *Aj_ we denote the dual vector in the x-y plane, that is in 
a Cartesian basis *A^ — A^, — —A^. The free Maxwell equations are then 

v_L • i?_L + = 0, v_L • + a^i; = , (7.1) 

V_L-*B_L-5t£; = 0, V_L +9*5 = 0, (7.2) 

d,B^~dt*E^=V^B, d,E^+dt*B^+ = V^E. (7.3) 

Locally, of course, the wave equation follows for each of the components, in par- 
ticular 

[a2-af-Aji? = o, [a2_af-Aji? = o. (7.4) 

Global extension to the whole region demands some boundary conditions for A_l 
in each case. We assume that the Dirichlet extension for the electric case and the 
Neumann extension for the magnetic case have been chosen. Assume moreover 
that /j^ Bdxdy = 0, that is B is orthogonal to constants in L^(Sl). Then one can 
represent these fields as 

B = -Ai*™, S=-A_L*e, (7.5) 

where Vfm is assumed to be orthogonal to constants, and the fields '^m and are 
thus uniquely determined by B and E. Equations (|7.4|) can be now expressed as 

[d? - dl ~ A^]vI/„. = , \dl - dl - Ai]*e = , (7.6) 

with appropriate boundary conditions. It is now easy to see that setting 

By_ = Vi9,*™ + *Viat*e , Ey_ = V^a.^e - *Vi9t*„ (7.7) 

one solves the complete set of Maxwell equations. The transversal fields satisfy 
boundary conditions 

n-i3_L = 0, *n-£;j_=0, (7.8) 
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where n is a vector in the x-y plane, orthogonal to dfl. Let us now add to B the 
omitted part B' independent of x± , and to B± and E± additional fields B'± and 
E' ± respectively in order to investigate uniqueness. We want to keep boundary 
conditions (|7.8|l . One finds then by integrating the first equation in H7.1|l and the 
second equation in H7.2|) over f2 that B' is in fact a constant, which we exclude. 
Then it is easy to show that the only solution for the remaining fields B'± and 
E' j_ is zero. 

Our boundary conditions are those of a perfect conductor, but this is not nec- 
essary in all its details. For instance, one could interchange the roles of the electric 
and magnetic fields. More generally, the boundary conditions should eliminate 
i_L-independent fields. However, for definiteness we keep our choice. 

The electromagnetic potential {(p, A±,A) producing fields ()7.5|l . H7.7|l may be 
chosen as 

if = -9,*e , Al_ = *V±*,„ , A = 9t*e . (7.9) 

However, the model is defined in terms of gauge-independent quantities and the 
potential is a purely auxiliary field. 

Using the boundary conditions one finds the total energy of the field in terms 
of and 



- [ [E^ + B^](t,x)d^x 
8vr J 

= 8^ E f[dt'^s{-^±)dt^s + d,^s{-A^)d,^s + (-A^^sf]{t,x)d''x. 

s—e^m 

(7.10) 

This brings us to the following formulation of the classical dynamics of our 
system. We use for this formulation variables Vg and introduced below, which 
should be thought of as supplying Cauchy data for fields and respectively, 
i.e. Vs = {'^ s , dt"^ s)\t=vn s — e,m. The electromagnetic field is derived from these 
variables. 

Let L'^{il)± denote the subspace of i^(^) orthogonal to constants. 



hi = y — and = y —A^, where is the Dirichlet selfadjoint Laplacian 

on Lg(Sl), and A'^ is the Neumann selfadjoint Laplacian on L^{n)±. Generally, 
let C be a selfadjoint operator on a Hilbert space H. By using a spectral represen- 
tation of C it is then easy to see the following: if C is positive and has a bounded 
inverse, then He = T^{C) is a Hilbert space with respect to the scalar product 
(iy9, i!)c = {Cip, Cip). Moreover, the restriction of C to I?(C^) is a selfadjoint, pos- 
itive operator in Tic, with spectrum equal to that of G. We apply this statements 
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to the operators and hy^ defined above, and denote 

ni^vih^), i^,^)^, = ^{hf^,hf^jj), (7.11) 

7^T=P(/i^), iv,ij)^^^^ihfip,h^^) (7.12) 

7^^ = L^(]R), 7^'= = 7^l ® 7^^ , 7^'" = (g) , (7.13) 

and the scalar product in 7?.* is denoted by (. , .)s. Moreover, we denote by hj_ and 
/i™ the selfadjoint operators in 7?.^ and TZ^ determined in the above way by 
and respectively, and then set 

= ^{hl^id)^ + {id®h,y, /i" = y^(/i^®id)2 + (id®/i,)2, (7.14) 

where hz = \J —d^. The space of our model is now 

L = Cr-®Z\ =V{h')®n' , s^m,e. (7.15) 

We write Vs = Vs(BUs €i C , s — m, e. The spaces are equipped with symplectic 
forms 

a%V,,V:) = i<,^s)s-ivs,<)s. (7.16) 
The Hamiltonian of the system is then 

n{Vrn,Ve)=n"'{V^)+n''{Ve), W{Vs) ^]^[{Us,U,)s + {h'^V,,h'v,),] , (7.17) 

and the evolution in each of the spaces is independently determined by 
as a symplectic transformation, as discussed in Sec. T3. One can now show that 
for sufficiently regular fields this reproduces the evolution equations 1)7. 6|l and the 
total energy of the field (|7.1()|l . 

It is now evident that our classical system is described by the di- 
rect sum of two systems of the type discussed in Section El If we put 
= {—Lx,Lx) x{—Ly,Ly), then /i^ and /i™ have the spectrum of the type (D) 
and (N) respectively. The quantization of each independent part of the system 
follows the same lines as before. The algebra of quantum variables of the entire 
system is then the C* -tensor product of the two Weyl algebras. The modification of 
the dynamics brought about by the modified boundary conditions is implemented 
by the replacement of H7.14|l by 

K = ^{hi®idr + {id®hi,r, K^ = ^{h^®idY + {id®hf,)\ (7.18) 

where /ij^, s = e,m, are constructed as in (|5.1I) . with hf^ determined by Neumann 
conditions for the case s = e and by the Dirichlet conditions for the case s = m 
(which imitates perfect conductor conditions) . The representations needed for the 
discussion of the Casimir effect are constructed as tensor products of the two 
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subsystems, and then the energy observable is the sum of the energies of the 
subsystems. Thus, in particular, the Casimir energy is the sum of the Casimir 
energies of the two subsystems, which proves our claim formulated in Introduction. 

Our strategy in the above analysis was to formulate the classical dynamics in 
terms of unconstrained and gauge-independent variables, and then quantize. For 
our purposes we did not need the quantum version of the electromagnetic potential 
or field itself. However, for the completeness we briefly sketch their formulation. 
We obtain the electromagnetic field from the potential, and for the latter we keep 
the quantum version of the formulas (|7.9|l . Recall that the fields independent of 
x± are to be excluded. Moreover, we do not expect the formulas for potential 
(or electromagnetic field) to extend to the boundary of O. Thus we take for our 
test function space Aj^I^n, where denotes the space of infinitely differentiable 
functions with compact support contained in f2 x K. Denote by and u'} heuristic 
quantized versions of Vs and Us respectively. Then motivated by H7.9|l we put for 

AiA^f) = J u'iA^fd'x = -M^^l, /)e , (7.19) 

A{A^f) = - j vl*^^A^fd'x^4n{vl,*^^fU. 

Now, the precise meaning of v'^ and is given by 

<^>s{Vs)^iv2,u,)s + iu'^,v,),, (7.20) 

where $s(K,) are quantum fields as described in Section 1-3. Therefore the precise 
meaning of H7.19|l is 

(^(Aj_/) = -4^$e(0,a./), A(Ai/) = -4^$e(/,0), (7.21) 
i*i(Ai/) =47r$„(0,*Vi/). (7.22) 

The electromagnetic field is then easily found: 

E{A^) = 47r<i>e(0, [A^ - (KJ^ + hl]f) , 
B(Ai/) =4^$„(0,Ai/), 

E^iA^f) = -47r$e(0, V^a,/) -4^$„(*Vi/,0) , 
B^iA^f) = 4^$e(*Vi/, 0) - 4^$„,(0, V^a,/) . 



(7.23) 



The formula for E depends on dynamics; for the free field h^^ should be replaced 
by hz, and then E{A±f) = $e(0, A_l/). Using the commutation relations 

[<i>.(y,),$,(K)] = *a^(T4,K)id (7.24) 
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one finds equal time commutators 



[B{A±f),E±iA±g)]=47n J iA±fyV±{A±g)d^x, 
[E{A^f),B^iA^g)] = -Ani f {A^fyV ^(A^g) d^'x 




(7.25) 



where e*-' is the antisymmetric symbol with e^^ = 1, and all other commutators 
vanish. For the free field the last term in the third relation vanishes and the commu- 
tators reproduce the usual quantization scheme. Also, the free Maxwell equations 
then hold for fields ifT?^ . In the presence of the modified plates the commutators 
cannot remain unchanged, as this would violate the constraints; the last term in 
the third commutator above takes care of that. The modification is nonlocal - this 
reflects the nonlocality of our model {{h%^Y is nonlocal). The pair of homogeneous 
Maxwell equations is still satisfied for fields (|7.23|l ~ by definition, which may be 
regarded as gauge invariance - while the r.h. sides of the other two give the sources. 
We note, moreover, that the identification of the full electromagnetic field in the 
interacting case is subject to some arbitrariness. The sources are linear functionals 
of the same fields, so one could shift some part of the electromagnetic fields to 
the r.h. sides of inhomogeneous equations to contribute to sources (but provided 
the homogeneous equations are conserved under the operation) . Our choice of the 
interpretation is the most simple one. 

We end with a summary of the physical meaning of our results. We believe 
that our analysis places the Casimir effect firmly within standard quantum theory. 
On the other hand it also shows why less conscious traditional quantum field 
formulations suffer from difficulties. The models analysed here are well defined and 
are free from the usual anomalies. The Casimir energy can be rigorously calculated 
for them as the expectation value of the positive free field energy operator in 
the ground state enforced by the environment (in accordance with the analysis 
presented in P). 

In the case of the electromagnetic field our results may have direct physical 
application. It follows from the discussion following (Asym) in Section that one 
should not expect a universal law for the total Casimir force between modified 
parallel plates. The original Casimir expression — 7r^/720a^ constitutes the third 
order term of the expansion of the energy in inverse powers of separation a, and 
cannot give a dominant contribution to the force for all values of this parameter. 
In particular, energy is positive, so the term is dominated by other nonconstant 
contributions for small a. On the other end of the range of a we have found that 
in our models typical fall-of of Casimir energy is governed by a term of order 1/a^. 
This prediction may be of more special character, but we believe that if it is indeed, 
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then disappearance of such a term in any other models needs further explanation. 
What we find in our setting is that for a wide range of models the coefficient at 
this term is positive, so the force at sufficiently large a becomes repulsive. 

Eventually, these predictions must be checked against experiment. However, 
although we have witnessed impressive progress in precision measurements of the 
Casimir force in recent years, and the existence of the force is now beyond doubt, 
its detailed form seems to us less certain. Due to experimental difficulties a mea- 
surement of the force between parallel plates has been reported only relatively 
recently p]. The results seem to confirm the original Casimir formula. However, 
we think it is to early to accept this as "the whole truth" on this force. First of all, 
the measurements span a limited range of the distance parameter a. As explained 
above we have very strong reasons to believe that the Casimir formula cannot hold 
for a tending to zero. The behaviour of the force for large a could also bring devi- 
ations. Second, the data analysis seems to be oriented at testing the coefficient at 
the 1/a* term of the measured force. To a large extent it involves peeling off other 
influences, while keeping the ~ l/a"* formula for the backreaction force; the best 
fit then confirms the original Casimir coefficient at this term. When seen from that 
angle the result does not contradict our predictions. We think that more extensive 
experimental tests of the functional form of the force are needed. 



Appendices 

A Sums and integrals 

In this appendix we gather a handful of approximation formulas connecting sums 
with integrals. Some of these results are well-known, but we give them the form 
needed in the main text. In particular, the Euler-Maclaurin expansion is usually 
formulated for analytical functions only, with no estimates on the rest. 

(i)A Let / : (0, oo) i— > {0,oo) be a non-increasing, continuous function. Then 

lim /(fee) = / f{u)du, (A.l) 

'^^ k=l 

/>oo 

hm e^ey J] f{khl + l^D = - / f{u) du , (A.2) 
<^x,e,-0 ^^^^^ 4 7o 

where the equalities hold both for finite and infinite integral on the right 
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hand side. If in addition 
»i 

f{u) du < oo then lim uf(u) — , (A. 3) 







f{u)du<oo then Hmu/(u) = 0, (A. 4) 

u — >oo 

poo ^ 

/ f{u)du<oo then lim ^ /(fc^e^) = q . (A.5) 



fc=l 



(ii)A Euler-Maclaurin expansion 

Let / be a complex function in C"((0, oo)) for some n G N and such that 



< C30, (A.6) 
lim/(™)(u)==0 for m = 0, l,3,5,...,<n-2. (A.7) 



Then the following identity holds 
lirn^ |/(0) + e J2 like) - / f{u) du 



(A.; 

n-l ^ 



where 



ml 

■m=2 



1 /■ 

(A.IO) 



6„(fc + s) - 6„Gs) =B„(1 - s) - iB„(0) - iB„(l) 



for s e (0, 1) , fc G N , 

%n,e|<C„/ (A.ll) 







Here are Bernoulli polynomials, i?„ = S„(0) are Bernoulli con- 

stants and c„ = ^i^-x l-Bri(s) — S„|. 

se(o,i> 

/"OO 

In particular, if in addition the integral / f{u) du converges then also the 

Jo 

oo 

sum /(fee) does (neither needs to converge absolutely), and in that case 

k=l 

the identity can be written as 

oo „oo n— 1 j-^ 

f{0) + eJ2fike)^ /(ii)d7.- ^e™^/("-i)(0) (A.12) 

fe=l •'^ m=2 
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A.l Proof of (i)A 

The first identity (|A.1I) follows immediately from Riemann approximations for a 
positive, non-increasing, continuous function 

1*00 />oo 

/ f{u)du<eY,fike)< f{u)du. (A.13) 
Je Jo 

The application of these inequalities to /(a^ + Z?^) as a function of a and /3 leads 
to 

/ f{a' + p') da df3 < e^ey ^ /(fc^e^ ^ ;2^2) 

/•OO /'OO pOO 

< / f{a^ + (3^)dadf3=j f{u)du, (A.14) 

which proves the identity ljA.211 . 

The function / is non-increasing, so uf[u) < 2 / f{v)dv, which implies 

Ju/2 

the properties (|A.3|I and l|A.4|l . By the first of these limits the first term in the 
sum in Eq. ljA.5|l tends to zero. For the rest of the sum, for e < 1, we have 



k=2 ^ •''^ 



f{n) , 
— ^ du 



<-Jj{u)du+^j f{u)du, (A.15) 

where the first inequality follows by the Riemann approximation of the integral. 
This is sufficient to conclude that Eq. (|A.5|I holds. 

A.2 Proof of (ii)A 

One first notes that for n G N the following identity holds 

N-l , 



fc=l 

n-l 

e™^[/(™-i)(0) - /("-i)(^e)] (A.16) 

m=2 
n N-l 

-^E / /^"nfce + s)[5n(l-(sA))-5^"(0)-^Bn(l)]rfs 
fc=o 



This is shown by induction with respect to n, with the use of integration 
by parts in the integrals in the second sum on the r.h. side (remember that 
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■j^Bn{u) = nBn-i{u)). The change of integration variable u = fee + s in each 
of the integrals in this sum puts the second line on the r.h. side into the form 

r / f ''Hu)bniu/e)du. 

nl Jq 

Functions &„ are measurable and bounded, so the estimate (|A.lip follows. There- 
fore, if the assumptions are satisfied then the r.h. side of (|A.16|1 . and the term 
{e/2)f{Ne) on the l.h. side converge for N oo. All the statements of the thesis 
are now readily seen. 

B Hilbert- Schmidt properties of operators Aa{u) 

In this appendix and hza are arbitrary selfadjoint, positive operators with 
domains ^{hz) and 'D{hza) respectively, and operators h{u) and ha{u) are defined 
by Eqs. H2.6|l and H2.7|l respectively (we keep this notation to make application 
of the results in the main text obvious, but we do not need any further general 
restrictions). Let / and g be real (or complex) measurable functions on (0,cx)) 
such that /(u), g{u)^ uf{u) and ug{u) are all bounded on the domain. Then the 
expressions 

f{hza)^V{u)g{hz) , where formally 

AP'^iu)^h-P{u)Aaiu)h-'^{u), Aaiu) ^ ha{u) - h{u) , 

define bounded operators for all p,q > and all u £ (0, oo) {u e (0, oo) if 
p = q = 0), although in general expressions A^'''(u) need not make sense by 
themselves. We denote by Pb and PaB the spectral projectors determined by 
and hza respectively, projecting onto the set B indicated in the subscripts. In par- 
ticular, the intervals in the symbols P{Pi,p2) ^^'^ -^0(71 ,72) below are chosen inside 
(0, 00) (positivity of operators). Symbol || A||hs will be used for the Hilbert-Schmidt 
norm of a HS operator A, i.e. ||A||hs = v^Trj^dM]. 

We shall show in this appendix that under these assumptions the following 
identities and implications hold. The principal results, which are needed in the 
main text, are contained in (iv)B and (v)b below. 

(i) B For each u > and v > the following integral representation holds in the 

uniform sense 

f{hza)Aa{u)g{hz) = - r [K{v)Aa{v) + Aa{v)h{v)\ i^Jt^dt . 

(B.2) 

(ii) B If /'a(-y^.^2)AP''?(u)P(^j^^2) is ^ HS operator for a given u = 1; > 0, then it is 

also a HS operator for u G (0, 00), and 
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)^a''(w)-f'(/3i,/32)llHS < ll^a(7i,72)^a''(«)^(/3i,/32)ll 



HS 



\Jll \J(il + ''^ I 1 ioTu>v, (B.3) 



(iii) B If 

Pa{^i,'i2)^'i (^)^a (^)-^(/3i ,/32) ^nd P^^^^ ^y^) ('-') (''^)^(/3i ,;92) are HS op- 
erators for a given u > 0, then Fa<7i,72)^a(w)-P(/3i,/32) is HS for u G (0,00), 
and 

X niax{||P„(^^,^^)/iy2(t;)A„(w)P(^^_^^)||„s,||P„(T,,,^,)A„(t;)/iV2(^)p^^^^^^^||^^|_ 

(B.4) 

(iv) B Let A2'*(w) be a HS operator for a given u = v > 0. Then it is also a HS 

operator for u e (0,oo), and 

IIA^NIIhs < ||Ar(t^)||Hs fovu>v, (B.5) 

II Ar WIIhs < (I) II Ar («)IIhs for < « , (B.6) 

hm||ArWI|Hs = 0. (B.7) 

The operator function (0, 00) 9 u 1— > A2'5(u) is continuously differentiable in 
the HS-norm sense and the formal differentiation yields the correct result. 
The operator AP'«(0) is HS iff ||AP'9(u)||hs < const for u G (0,oo) and then 
||Ar(0)||„s= lim||AS'''(u)||Hs. 

The theorem remains valid upon replacement of the operator A^J''(u) with 
PaB^a''^{u)Pc, for any measurable sets B,C. 

(v) b Let Aa(0) be a HS operator. If Aq('u)/i^/^(u) is HS for a given u = v >0, 

then all values of the operator fmictions (0, 00) 3 u Aa{u)h^^'^{u) and 
1/2 

(0,00) B u ha {u)Aa{u) are also HS operators. Both functions are con- 
tinuous in the HS-norm on their domain and continuously differentiable in 
HS-norm sense on (0,oo). Moreover, 

hm \\hl/\u)Aaiu)\Us = , lim \\A,{u)h'^\u)\\^s = ■ 



The theorem remains valid upon replacement of the operator Ao(u) with 
PaB^a{u)Pc, for any measurable sets B,C. 
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B.l Proof of (i)B 

For a positive real number a one has 1 — — 

TT 

spectral representation of hz one shows that for u > 



dt 



Using this in the 



h{u)g{hz) 



1 



{hi + u)g{hz) dt 



and the integral on the r.h. side converges uniformly (in norm). With a similar 
representation of f{hza)ha{u) we have a uniformly convergent representation 



f{hza)Aa{u)g{hz) = 3 / fihza) 

Using the formal relation 



hl„ +t hl+t 



dt 



\/t — u 



hl„ + 1 hl + t 



{t~u) 



[/la(w)AJt.) + Aa(w)/l(v)] 



hl+t' 



which becomes a correct identity when placed between f{hza) and g{hz), one 
arrives at Eq. (|B.2|) . 



B.2 Proof of (ii)B 

If one multiplies Eq. (|B.2|I by h~P{u) from the left and by h~'^{u) from the 
right {u > 0), one obtains a similar identity with AJJ'* replacing Aa on both 
sides, and the integrand on the r.h. side multiplied by hP{v)h~^{u) from the 
left and by h'^{v)h~'^{u) from the right. We put f{hza) = ^'0(71 ,72) (^za) and 
5(^2) — Pn3ijh)i^z) in this identity, and estimate the HS norm of the l.h. side. 
We find 



\\Pa (71,72 

)^a'^i'^)P{l3ul32)\\iiS < \\Pa{ 71,72 

X (ll^a(7i,72)'ia(w)|| + W P{f3ul32)Hv)\\) 

^ j_ Pa{li,l2)i^za) 
L + i 

where we have pulled the HS norm sign under the integral and used the fact 
that IIABCIIhs < ll^ll II-B||hs||C'||- The second line in this estimate is bounded by 
VtI+w + \A^2+^' fi^'st factor in the third line by max{l, the 
second factor in the third line by max{l, (w/u)*/^}, and the fourth line by 

1 ^Jt^dt _ 1 

^ Ju hT+t)W+^^^W+^+VW+^' 

which ends the proof of (|B.3|I . 



P 



){hz) 



hl+t 



\/t — udt , 
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B.3 Proof of (m)B 

We put f{h^a) = Pa{it,i2){hza) and g{h^) = P{i3i,fi2){hz) in the identity (|ij.2|). 
Using the method of the last proof we find 

ll-Pa(7i,72)^a(")^'(/3i,/52)l|HS 

— (11^0(71,72) ^y^(")lll|-Pa(7i,72)^a''^(«)^a(w)P(;3i,/32)llHS 

./92)IIhs||^'(/3i,/32)^^^^(^')II) 



-^0(71,72) i^z 



which leads to the estimate (IB.4I 



hl+t 



y/t — udt , 



B.4 Proof of (iv)B 

We set in the estimate (IB.3II 



where fc, Z = 0, 1, . . ., and denote 



au {t,v,u) 



\/kT + T + V + \JIt + r + ' 
^/kr + u + ^/It + u 



Also, we introduce Pa^k,r = Pa{^.Vk¥+^)^ ^'^^ = P{VT^,VVF+^y '^'^^'^ 
bound we have 



< |lPa,fc,rAP^'?(u)Pi,^||Hsafci(r,W,M) X |^^^(p+g)/2 



for u> V . 
for u < V . 



If V < u then we choose r such that v + t < u, and then aki{T, v, u) < 1. If u > u, 
then aki(T,v,u) < [{y + t) /uY^"^ . Moreover, there is limu^oo — 0- It is 

now sufficient to observe that 



2 

HS 



to be able to conclude that Eqs. (|B.5p and IB.Tp hold. For u < w one ob- 
tains |lAP'«(u)||2<, < ||AP'9(w)||23 [(w + r)/w](u/u)P+« for each r > 0, which 
leads to the estimate (IB.6|I . Also, it should be clear that the replacement 
AP'9(.) ^ PaBAP'«(.)Pc poses no difficulties. 
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To save space we prove the remaining statements only in the case p = q = 0, 
the generahzation to nonzero p, q and/or added projections PaB, Pc is then easily 
obtained. We note first that 



h{t) - h{u) 



t — u 



h{t) + h{u) 



and similarly for ha{-), are bounded operators. Thus if A(j(t), Ao(u) are bounded 
then 



A„(t)-A„(u) = -(i-u)- 



1 



ha{t) + ha{u) 

If on top of that Aa(f) and Aa(u) are HS, then 

\t-u\ 



(Aa(t) + A„(u)) 



1 



h{t) + h{u) 



|A„(i)- A„(u)||hs < 



(Vi + V^f 

In a similar way one shows now that for t,u> Q there is 



Aa(i)||„s + ||A„(u)||Hs). 



Aa(i)-A„(w) 1 ^1 1 



< const - 



|A„(t)||HS + ||A„(u)||Hs)- 



Therefore, if the assumption is satisfied (with w > 0), then the operator function 
(0,oo) 9 u ^ ^a{u) is HS-difFcrcntiablc, and its derivative is the HS-continuous 
operator function — ^h~^ {u)Aa{u)h~^ (u) . 

If Ao(0) is HS, then the function (0,00) B u 1-^ ||Aa(u)||Hs is decreasing, 
hence ||Aa(0)||Hs > lii^ l|Aa(i)llHs > ||Aa(u)||Hs for u G (0, 00). Conversely, if the 
decreasing function (0,oo) 9 u 1— > ||Aa(u)||Hs is bounded, then 

N 

00 > lim ||A„(t)||2^ > ||A„(^.)||L > E l|AaW<^n||^ 
where {ipn} is any orthonormal basis. But 



n=l 



|A„(w)-A„(0)|| 



ha{u) + hza h{u) + h^_ 



so taking the limit u \ 0, followed by ^ 00, we have lim ||Aa(t)||Hs > 



I A, 



. In conjunction with the opposite inequality obtained above this be- 



comes the desired equality. 
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B.5 Proof of (v)b 

If Ao(0) is HS, then by (iv)B also Aa{v) is HS. Thus if in addition Aa{v)h^/'^{v) 
is HS, then by observing the identity 

Aaiv)haiv)Aaiv) = AI{v) + Aaiv)h{v) Aa{v) 

we learn that hi^'^ (v)Aa{v) is also HS. The proofs of the theorem for Aa{u)h^^^{u) 
and ha {u)Aa{u) are similar, we take the first of these functions. Using the esti- 
mate (IB.4|) we have 

\\Pa,k,rAaiu)h'/^{u)Pi^,\\l, < (Zt + T + u) || A J^,)p,^^ || 2 ^ 

<cli{T,V,u)(^\\Pa,kAl^Hv)Aa{v)Pi,r\\l, + \\Pa,k^rAa{v)^^^^ 

where 

, ,w4(fcr + T + w)i/4 + (;r + T + i;)i/4 

CkliT, V, U) = {IT + T + V 

\/ KT + U + \/ IT + U 

< o|- '^ + '" \i/4 (fcr + T + w)!/'' + {It + t + 



2 

HS y ' 



' (fcr + u)i/4 + (/r + u)i/4 



Thus Aa{u)h^^'^{u) is a HS operator for all u e (0, oo). Also, lim Cki{T,v,u) = 

u — >oo 

and Cfe;(T, V, u) < 2^1^ for u > r + w, so || Aa(w)/i-'-/^(it)||Hs ^ for u oo. 

The HS-differentiability on (0, cxd) is proved similarly as in (iv)B. It remains 
to investigate Aa(0)/i^/^(0). Using the methods of the previous proof we find 

ha(u) + ha[0) h(u) + h{0) 

+ A„(0) 



[h^^iu) + h^/^{0)][h{u) + h{0)] 



As the r.h. side and the first term on the l.h. side are HS operators, the operator 
Aa(0)/ii/2(0) is also HS. Estimating the HS -norm of the r.h. side we find 

\\Aaiu)h^/'iu) - Aam^/HO)\U < ii^/'*(2||A„(0)||Hs + ||A„.(u)||„s) , 

which proves the missing HS-continuity at zero. Finally, it is now easy to convince 
oneself that the replacement of Aa(.) with PaBAa{.)Pc poses no difficulties. 
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C Definitions, estimates and identities 

In this appendix we introduce some denotations and prove some results needed for 
the evaluation of the integral kernel of the operator hza — as defined by (|5.1|l . 

Let be a real (or complex), bounded function in C^'^^((0, +00)) for some 
N >0, such that |F(^+2)(p)| < const {p + i)-(^+2+7) for some 7 E (0, 1). Then 
also 

\F^"'\p)\ < const {p + 1 < 71 < iV + 2 , 

\F{p)-Fo,\< const {p+l)-\ ^^'^^ 

where i^oo is a constant (limit value at infinity). We define the following functions. 
For {q,p) e R2 \ {0, 0} we denote 

Fi\q\)~F{\p\) 

>^F[q,P) = -, w , X , (C.2) 

(q-pjiq+P) 

so that 

{q,p) ^ Af (p, g) = Af (-<7, p) = Af {q,-p)- (C.3) 

We denote for p e M 

/ + OO POO 
XF{q,p)dq^2p Xpiq^p) dq = -Af{~p) , (C.4) 
-00 Jo 

+ 00 C30 

^fAp) = -P^ X! •^^'(^^'P) = P'^^f{0:P) + 2pe^Xp(ke,p) = -Ap^^{-p) , 

k— — oo k—1 

(C.5) 

AO,,,(p)=peAF(0,p) = e^^M^:ZM, A>,(p) = 2pe V AF(fce,p) , (C.6) 
and for p ^ p': 

, Af{p)~Af{p') ^ . ,x 

Xf(p,P ) = -. = Xf(P ,P) = Xf{-P, -P ) , (C.7) 

P-P 

, ^fAp)~^fAp') / / n / /X ,r'Q\ 

xfAp^p ) = } = xfAp iP) = XF.A-p, -p ) ■ (^-■•S) 

p-p 

, „ A°,^,(p)-A"^,(y) ^ A^,(p)-A^,,(pO 

XF.eP^p = r- ' Xf,Ap^p) = — } • (C-9) 

p — p p — p 

With these assumptions and denotations we have the following results. 

(i)c The function Af is in C^+HlR+), aU derivatives A^"'"'' for < m + n < iV+1 
extend to continuous functions on (0, 00)^ \ {0, 0} and satisfy the estimates 
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for g +p > 1 : |A^'"^(g,p)| < const < 



1 

, 9 > P 



|^(g + l)m+7(p+l)n+2 ' ^"P' 



The integrals A^'°''(g,p) are in C^+i-™(IR+), the differentiation may 
be carried out under the integral sign, and one has the estimates: 

roo fdlogpl + l), m + n = 0, 

forp<l: / \xf'''\q,p)\dq< const I 1 ^ (C.12) 

Jo y-;^^ m + n>l, 

1*00 

for p > 1 : / |A^'"''(g,p)| < const 
Jo 



1 

m = , 



,n+l+7 ' 

-TTT , m > 1 



PV-' (C.13) 



(ii)c The functions Aj?, AJ^,, and hPp^^ are in C^+^(M+) and satisfy the estimates 



{p{\\ogp\ + l), n = 0, 

|A(,")(p)| , |A;,5)(p)| < const I (| logp| + 1) , n = 1 , 
forp<l: (C.14) 



|A^(:^(p)| < const e, 



for p > 1 : \K^-\p)\ , |A;5)(p)| < ^ , |AO,,,(p)| < ^ . (C.15) 

The functions xf and XF,e are in C^((R \ {0})^) and in this domain satisfy 
the estimates 

for < IpI < 1, < \p'\ < 1 : 

\xf{p,p')\, |xF,.fep')l <const(|logH| + |log|p'|| + l), (C.16) 

for < IpI < 1 , < b'l < 1 , pp' > : \x%,eiP,P')\ < const e , (C.17) 

const c 

for < bl < 1 , < b'l < 1 , pp' < : |x^,,(p,p')| < ' C^'^^) 



for (p,p')GR'\(-l,l>': 



\xf{p,p')\ , \XF,eiP^P')\ < const < 



, |P|<|P'I 



+ m(b'| + i)' (C.19) 
|p|>|p'l- 
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Moreover, for \p\, \p'\ > 1 there is 



\x^r''\p,p')\, |xte'"'(p,p')l< const ^ 



(m,n) / 



n+1 



(|p| + l)m+7^(|p/| + l) 

l(bl + i)™+i(b'| + i)"+T 



bl>b'l- 

(C.20) 



(iii)c The following identity is satisfied in the distributional sense: 



/+00 
-oo 



(p — q ± iO){p' — q ± lOj p — P 

On the l.h. side the distribution [{p — q± iO){p' — q± is first applied 

to a smooth function of compact support f{p,p'), and the result integrated 
as indicated; the r.h. side is multiplied by f{p,p') and integrated over dp dp'. 
Note that in the integral on the l.h. side the signs in front of iO must match. 

For p, p' ^ ke there is 



+00 

^ E 

fe=— 00 



+00 

^ E 

k=—oo 



F{\k\e) 



(p — ke)(p' — ke) 



(C.22) 



= xfAp^p') 
(-i)^nifck) 

{p — ke) {p' — ke) 



p-p 



-[F{\p\) cot{wple) - F{\p'\) cot(7rp7e)] , 



XF,2e(P,p') - XF,e{P,p') 



P-P 



F{\P\) 



n\p'\) 



(C.23) 



sin(7rp/e) sin(7rp'/e) 



In consequence further identities follow for all p, p' G M: 



+•00 



sin[6(p — ke)] sin[&(p' — ke)] 



= COs[b{p - p')]XfMp') + COs[b{p + p')]\xF,e{P^p') - XF,2e{P^p')] ^^"^^^ 



ip - p') 

^' ' '^^ ' {p + ke){p' - ke) 



k=—oo 



■ cos[b{p - p')]xfAp^ -P') - cos[6(p + p')] [xF,e(p, -p') - XfMP' -p')] 

TT sin[6(p — p') 



{p + P') 



m\p\)-Fi\p'\], 



(C.25) 
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C.l Proof of (i)c 

For q,p > one has 

F{q) - Fip) 



q-p 



F^^\qt + p(l ~ t)) dt , 



SO by (|C.1|I the differential properties of A^'"^ are satisfied. For p + q < 1 the 



(C.26) 



derivatives of (|C.26|) arc bounded, so the estimate (|C.10p is also true. To prove 
the estimate (jC.llll it is sufficient to assume that q+ p>\ and q> p (due to the 
symmetry oi Xp), which implies q + p> {1 + q)/2 and 



1 



< 



const 



(C.27) 



We consider two cases q < 3p and q > 3p separately. In the first of these regions 
one has p > q/3 > {q+p)/6> {q + 1)/12 and 



grgs 



F{q) - F{p) 



< 



q-p 

< const 



r+s+l) 



{qt+p{l~t))\dt 



dt 



< 



const 



(C.28) 



ip + {q - p)tY+'+^+^ - {q + lY+'+^+-f 



so |A*^™'"''((7,p)| < const {q + i)-(™+"+2+7) ^jjig region, which complies with the 
estimate (|C.llll . In the second region one has q — p = {q+p + q — 3p)/2 > {q+l)/4:, 
so taking into account Eq. 1C.27|I one has in that region: 



d'^d'' 

1 P^q.p)^q+p) 



< 



const 



(C.29) 



Taking also into account the bounds (|C.ip . and similar ones with the argument p 
replaced by q, one confirms the estimate (jC.ll|l . 

The statements following estimates (ICllOp and are their simple con- 

sequences. The next two estimates (|C.12ll and (ICLSH are obtained from the pre- 
ceding two by elementary integration. 



C.2 Proof of (ii)c 

For Ap the differentiability properties and the estimates (jC.14|) and (jC.15|l are 
simple applications of the properties stated in the last sentence of (i)c- To prove 
the same facts for Ai?^e one has to estimate the sums by integrals, which can be 
done rather easily by observing that for each p > the r.h. sides of the bounds 
((Cno)l and fCnTjl can be glued together into a continuous decreasing function of q. 
For A^^ the estimates (IC.14II and (IC.15|) follow from (|C.26|) and HC.11|) respec- 
tively. The estimates (|C.20p in the enlarged region > 1/2 follow from the 
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bounds ljC.15|l (which stay valid for \p\ > 1/2) with the use of the method similar 
to that applied in the proof of (i)c- The estimates (|C.19|) are then easily extended 
to the whole region of application with the use of bounds HC.14() . HC.15|I (with 
n = 0). The remaining estimates also follow from these bounds by cutting the 
square (—1, 1)^ into four squares according to the signs of p and p' . 

C.3 Proof of (iii)c 

To prove ljC.21|l we first note that for complex z, z' with 52, Qz' 7^ there is 



z 



— 00 



(g-z)(9 + z) {q-z'){q + z') J _^ {q - z){q - z') 



(use evenness of the integrand on the r.h. side). Let 1 > a > a' > 0. Setting 
z — p ± ia, z' — p' ± ia' , denoting 



Fi\q\) - F{\p\) 
(q — {p± ia)){q +p±ia) 



and using the identity 

'■+°° dq 



±- 



{q — {p ± ia)){q + p ± ia) p±ia 



we get 



Fi\q\)dq A±"(p)-A±«'(p') Fi\p\) - Fi\p'\) 

— ± ITT- 



) {p — q i ia){p' — q ± ia') p — p'±i{a — a') p — p' ± i{a — a') 

(C.30) 

We multiply this identity by a Schwartz function f{p,p'), integrate, and take 
successive limits a' — > 0, a — > 0. Now, for the l.h. side we observe that the order 
of integration may be changed and the result represented as an integral over q of 
F(|(7|) multiphed by 



/- '^^w ^ ^ rr dp dp' . 
(p — q ± ia) (p' ~ q ± ia') 



By the well known techniques one shows that the consecutive limit exists and, 
moreover, the above expression is bounded by const {q^ + 1)~^. This is sufhcient 
to conclude that l.h. side of Eq. HC.21() is obtained in the limit. For the r.h. side 
of ljC.30|) one easily notes that its consecutive point-wise limit yields for p ^ p' 
the function on the r.h. side of ljC.21|) . Thus to complete the proof one only needs 
to show that in the limiting process the r.h. side of (|C.30p stays bounded by a 
function defining a distribution. This is immediate for the second term, as 



Fi\p\)~F{\p'\) 



p — p' ± i[a — a') 



< \p + p\ \Xf{p,p) 
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Also, the property is rather obvious for the Hmiting process a' in the first 
term (note that \Ap"{p)\ < {\p\ + 1) J \XF{q,p)\dq). Thus we are left with the 
function 

p — p' ± ia p — p' ± ia p — p' ± ia 

The second term on the r.h. side of this equality is bounded by \xf{p,p')\^ which 
is sufficient for our purpose, while for the first one it is sufhcient to estimate the 
function 

Apip)~AFip) f [Fi\q\) - F{\p\)] dq , . 2p±ta 

T- " / 7 / I ■ ^^/ — ] — — — ±^7rp >^f{p,P) 

±ia J [q - [p ± ia)){q + p ± la) p±ia (C 31) 

, /o , • ^ /■ >^F{q,p) - Xf{p,p) , 
+ p{2p±ia) / / , . Nw , I ■ ^ g^g- 

J [q - [p ± ia)){q + p ± la) 

The first two terms on the r.h. side of this identity are bounded respectively by 
/ \\F{q,p)\dq and 2T:\p\p[p,p)\, which fulfills our demands. For the third term we 
assume that p > Q (the case p < Q needs only obvious modifications), represent the 
integral as twice the integral over (0, oo), and consider the integration sets (0,p+l) 
and (p + 1, oo) separately. The term containing integration over the second set is 
bounded by 



log(l + 2p) 

qe{p+r,<x>) ^'="'^'7 - (p + 1)1+7 

The remaining part is bounded by 



(p + 1) log(l + 2p)( |Af(p,p)| + max |A_f((3',p)| )< const 
^ oe(p+i,oo) / 



2p 



>^F{q,p) - Xf{p,p) 



p+i 

dq<2p I j-^ I ' |A^'^")(s,p)|ds dq, 



where r< = imii{q,p}, r> = maji{q,p}. For p > 1 there is |A^'°^(s,p)| < constp ^ 
(cf. (|CnT|l l. which gives a sufficient estimate in this region. For p < 1 and q < p+l 
one has |Af(s,p)| < const {s+p)^"^ (by (|C.10|I . which may be used in this extended 
region). Thus the r.h. side is bounded by const log[(2p+l)/p], which ends the proof 

of irriiii . 

To prove (|C.22|) we first note that for p,p' ^ ke, there is 



But 

^^E (fcT^ll^ = A.,(p)-^(H).cot(.,A), 
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where we used the fohowing identity 

— 



9^,9 ^ 7rcot(7rx) , (a; ^ Z) 



(cf. m, formul a 1.421 (3)). This ends the proof of IHSHl- 

Formula (|C.23|) follows easily from the preceding one if one observes that 

. ., , . „ ^ F(\k\2e) ^ F(\k\e) 

its l.h. side may be written as 2e ) ^ , ^ , , , — — — e > , , , , , -— . 

(p - k2e){p' - k2e) (p - ke){p' - fee) 

Finally, to prove ljC.24|l and ljC.25|l one notes that 

2 sin[6(p - fee)] sin[6(p' - fee)] = cos[6(p - p')] - (-l)*^' cos[6(|7 + p')] 
and uses !^l22l and (Rl23l). 



D Some identities for the logarithmic derivative 
of the Gamma function ifj 

Two particular textbook representations of the function ip are of importance for 
us: 

1 r ^11 

^(z) = — + hm logTV - ^ — - , (D.l) 

1 r°° 

Tp{z) - log z + - = / v{s)e-^' ds, ^z>0, (D.2) 

where 







vis)^--^. (D.3) 

s — 1 

The function v is analytical in the complex plane outside the points z ~ 2fc7ri, 
k G Z \ {0}, where it has poles with principal values —1. Moreover, 



,W(0)._|^, |,W(,)|<__£i_, .eM, fc = 0,l,.... (D.4) 



Using these properties in the representation (|D.2p one finds by induction for m G N 
and ^z > the expansion: 

1 ™ f> 1 1 /"OO 

^W-l°g- = -^-Elf^ + ;5;^/ -''"H^)e-ds (D.5) 
k=i 

(remember that i?2fe+i — for fc G N). We denote 

w,^iz) = z'-L{z) _ logz + 1 + . (D.6) 

k=l ■' 
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Then for A > 0, m = 1, 2, . . ., one has the identity 
In particular, 

and by the apphcation of the operator Xdx: 

To prove Eq. (|d1^ one notes that z;(2™)(_s) ^ _^;(2m)(s) and uses Eq. ||D^ 
to obtain for A > 

io cos(2.A0-.„(t)d. = -y^ 7T72^^^- (7T72^X)^^^- 

(D-10) 

The proof is now completed for A > by integration in the complex s-plane along 
the rectangular contour with vertices at points ±^ + (2fc + l)ni followed by 
the successive limits ^ oo and k ^ oo. Finally, one uses the continuity in A of 
both sides of 1D.7II to extend the formula to A = 0. 
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